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' Abstract 
> ■ 

• We develop a gauge invariant canonical perturbation scheme for perturbations around symmetry 

I reduced sectors in generally covariant theories, such as general relativity. The central objects of 

. investigation are gauge invariant observables which encode the dynamics of the system. 
' We apply this scheme to perturbations around a homogeneous and isotropic sector (cosmology) 

' of general relativity. The background variables of this homogeneous and isotropic sector are treated 
fully dynamically which allows us to approximate the observables to arbitrary high order in a self- 

' consistent and fully gauge invariant manner. Methods to compute these observables are given. 
^ , The question of backreaction effects of inhomogeneities onto a homogeneous and isotropic ha,c]i- 

O^' ground can be addressed in this framework. We illustrate the latter by considering homogeneous but 

I ' anisotropic Bianchi-I cosmologies as perturbations around a homogeneous and isotropic sector. 
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1 Introduction 



\^ ' General relativity has two very challenging features: firstly the dynamics of the theory is highly non- 



linear, secondly general relativity is a diffeomorphism invariant and background independent theory. 
These two features make it very difficult to construct gauge invariant observables, that is to extract phys- 
ical predictions. Diffeomorphism invariance of the theory includes invariance under time reparametriza- 
tions, therefore observables have to be constants of motions. Hence finding gauge invariant observables is 
intimately related to solving the dynamics of the theory. But because of the highly non-linear structure 
of the theory it is quite hopeless to solve general relativity exactly. Indeed so far there are almosi0 no 
gauge invariant observables known [T]. 

Therefore we think that it is important to develop perturbation schemes in order to attack the dy- 
namics of the theory. The difficulty for a perturbation scheme in general relativity is again the control 
over the gauge dependence of the results , which becomes already quite challinging at second order (see 
for instance 

In this work we will develop a manifestly gauge invariant perturbation theory in the canonical frame- 
work. This perturbation theory allows the calculation of observables to an arbitrary high order. In an 
earlier work we introduced a calculation scheme for gauge invariant observables for perturbations around 
a fixed background spacetime, which in the canonical framework corresponds to perturbations around a 
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fixed phase space point. In this work we will rather consider an expansion around a whole symmetry re- 
duced sector of the theory. This sector will be treated non-perturbatively and could for instance describe 
homogeneous and isotropic fields (i.e. cosmology), spherically symmetric fields (i.e. black holes) or even 
midi-superspaces such as cylindrically symmetric gravitational waves. The latter example has infinitely 
many degrees of freedom, but is solvable [S] l4] . Note that perturbations around a fixed phase space point 
arise as a special case of the scheme developed here. 

We think that the methods developed here are useful for a quantum theory of gravity as well as for 
approaching classical dynamics, for instance for cosmological perturbation theory. 

In a quantum theory of gravity gauge invariant (that is physical) observables are the central object. 
The construction and interpretation of the observables of the theory is therefore one of the key open 
issues, which is also related to the problem of time, see for instance [3 El El O [H [10]. Approximation 
methods for gauge invariant observables provide an explicit way to construct observables and may also 
help to test (quantum) interpretation of these as well as to discuss phenomenological implications [3 [9] . 
The approximation scheme developed in this work shows explictly how to express local observables in a 
gauge invariant way and therefore indicates how one could understand quantum field theory on curved 
space time as an approximation to the full theory of quantum gravity. 

Using an approximation scheme for observables around a whole (symmetry reduced) sector of the 
theory allows one to explore properties of gauge invariant observables better than in a perturbation scheme 
around a fixed phase space point. First of all one can now incorporate results from symmetry reduced 
(exactly solvable) models. The degrees of freedom describing these sectors are treated non-perturbatively. 
Secondly, as we will see, this approach provides a gauge invariant (perturbative) description of dynamical 
processes also in closed universes, where for instance the notion of scattering amplitudes is problematic. 

Moreover since the early days of quantum gravity mini- and midi-superspace models gained by sym- 
metry reduction of the full theory played an important role, especially for describing quantum cosmology 
[12] and references therein). However because of the non-linear dynamics of the theory, which leads 
to a coupling of non-symmetric and symmetric modes, the question arises whether the results gained 
from these models are reliable |13j . For different attempts to derive symmetry-reduced quantized models 
from the quantized full theory, see for instance [HI [151 [16] . Since the coupling between non-symmetric 
and symmetric modes arises through the dynamics, we expect that the issues raised here become im- 
portant if one wants to match solutions to the Hamiltonian constraint in the symmetry reduced model 
with solutions to the Hamiltonian constraint in the full theory. We provide such a matching for classical 
gauge invariant observables, however the observables in the full theory have higher order corrections due 
to the non-symmetric modes. These corrections can be given explicitly, thus allowing for an estimation 
of errors, which arise if one ingnores the non-symmetric modes. 

Recently there is also growing interest to incorporate (linear) perturbations [17' into the framework of 
loop quantum cosmology |,18| . Here we hope that a clear gauge invariant formulation of perturbations of 
symmetry reduced models might help to understand for instance issues related to the gauge dependence 
of these approaches. Since the scheme developed in this work is consistent to any order, it proves that 
linear perturbations can be understood as lowest order of fully gauge invariant perturbations governed by 
non-linear dynamics. Similar gauge issues are present in recent developments in background independent 
approaches to the graviton propagator [19] . 

A gauge invariant canonical perturbation theory could be also fruitful in classical applications, such 
as second (and higher) order perturbation theory around cosmological solutions [20] or black holes. The 
main difficulty here is to control the gauge dependence of the results. This gauge dependence can be 
understood from the fact, that one has to identify spacetime points in the "physical" (non-symmetric) 
universe with spacetime points in the "background" universe, around which the perturbation is taken. 
This identification can be related to a choice of coordinates for the "physical" universe. One might wonder 
why we attempt to develop a perturbation theory in the canonical formalism, where one would expect 
the problem to be even worse due to the foliation for the "physical" and "background" universe one has 
to choose in the canonical framework. 

The resolution is that we use observables as central objects, i.e. we attempt to approximate directly a 
gauge invariant observable of the full theory and do not consider (the difference of) fields on two different 
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manifolds representing the perturbed and unperturbed spacetime. Observables in the canonical formalism 
correspond to phase space functions, gauge invariant observables are invariant under the action of the 
constraints (the gauge generators). 

The phase space of general relativity is just a representation of the space of all spacetimes (i.e. 
solutions of the Einstein equations). Gauge invariant phase space functions give the same value on 
spacetimes which are related by a diffeomorphism. Hence by considering gauge invariant phase space 
functions we do not need to worry about the identification process between points in the perturbed and 
unperturbed spacetime. 

In order to approximate gauge invariant phase space functions we have to declare which variables 
are to be considered small. This choice is done in such a way that the approximate observables coincide 
with the exact observables if evaluated on the symmetry reduced sector of the phase space. Indeed the 
zeroth order variables can be defined by an averaging procedure. First oder phase space functions vanish 
on symmetric spacetimes, higher order phase space functions are products of first order phase space 
functions. Note that the splitting of phase space variables into zeroth and first order is done on the gauge 
variant level. Generically a gauge invariant phase space function is a sum of terms of different order. 

This approach has similarities to the (1 + 3) covariant perturbation theory [21j which introduces a 
preferred timelike observer congruence. In this work we rather have to declare a part of the degrees of 
freedom, as for instance scalar fields or the longitudinal modes of the gravitational field, to be used as 
clocks, which might be a more general procedure. However we think that it should be possible to recover 
the (1 + 3) covariant perturbation theory if one manages to match the choice of clocks to the preferred 
observer congruence. On the general relation between observables in the canonical theory and "covariant" 
observables, see [7]. 

A key feature of this work is, that we keep the zeroth order variables as fully dynamical phase space 
variables and not just as parameters describing the background universe as one would do in a perturbation 
around a fixed phase space point. Indeed we have to keep the zeroth order variables as canonical variables 
to allow for a consistent gauge invariant framework to higher than linear order. Moreover this provides 
a very natural description for backreaction effects: these arise as higher order corrections to observables 
arising through averaging of (time evolved) phase space variables. Since this approach is gauge invariant 
it could shed some light on the discussion whether these backreactions are measurable effects or caused 
by a specific choice of gauge, see for instance HHHnilll]- As already mentioned we have to choose clocks, 
which define also the hypersurfaces (by physical criteria, e.g. by demanding that a scalar field is constant 
on these hypersurfaces) over which the averaging is performed. Therefore the observables describing the 
backreacion effects depend on the choice of clocks. However, as we will see, one can find relations between 
the gauge invariant observables corresponding to one choice of clocks and the gauge invariant observables 
corresponding to another choice of clocks. 

Let us shortly describe the main ideas of the approach developed here: We will approximate a special 
class of gauge invariant observables, known as complete observables [6l[25l[7]. The complete observables 
have the advantage that they describe the dynamics of the theory by giving the evolution of certain 
non-gauge invariant observables (the partial observables) with respect to other non-gauge invariant 
observables (the clock variables). In [25] methods to compute these complete observables were developed, 
in particular a power series was derived. 

This power series serves as a starting point for our approximation scheme. We will devide the phase 
space variables into two sets, namely variables of zeroth and variables of first order. This devision can be 
implemented by using a projection operator on the space of phase space functions which projects onto 
the symmetry reduced sector. First order variables are projected to zero. For cosmological applications 
this projection operator is defined by an averaging procedure. Now one can define approximate complete 
observables of order k by omitting in the power series all terms of order higher than k. 

Moreover one can define evolution equations for these complete observables and even gauge invariant 
functions, that generate this evolution (that is Hamiltonians) . As we will see the calculation of the 
complete observables up to a certain order can be cast into a form where "free propagation" (i.e. the 
linear propagation of the first order variables on the symmetry reduced background and the evolution 
of zeroth order variables in the symmetry reduced sector) , is perturbed by interaction processes between 
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first order fields and between first order fields and the zeroth order variables. 

In the next section [2] we will shortly summarize the necessary details concerning complete observables 
for a general gauge system (in a canonical description). Furthermore we will define gauge invariant 
Hamiltonians, which generate the physical time evolution with respect to the clocks. Section [3] introduces 
the approximation scheme for a general system around a symmetry reduced sector. There we also define 
approximate gauge invariant observables of a certain order and consider first properties of approximate 
complete observables. 

In section 2] we consider this scheme for perturbations around a sector describing isotropic and homo- 
geneous cosmologies. We derive and interpret the power series expressions for the complete observables. 
Then we consider the transformation between complete observables defined with respect to different 
choices of clocks (which can be understood as representing different families of observers) . Section [6] de- 
fines lapse and shift functions, which allows to compare this canonical approach to covariant approaches. 
Lapse and shift function determine foliations of spacetimes. These foliations are defined by the choice of 
clock variables, i.e. by physical conditions. 

Section [7] gives explicit examples for clock variables. We consider clocks built from the gravitational 
degrees of freedom and clocks defined by the matter field degrees of freedom. For the former case we can 
give a set of clock variables which is related to the longitudinal gauge (see for instance [17]), the dynamics 
of the corresponding complete observables is considered in appendix [b1 

In section |5] we consider the equations of motion for first order complete observables (associated to 
the scalar modes) and a method to find these complete observables. Then we consider the lowest order 
backreaction effects, that is second order complete observables associated to zeroth order functions. We 
evaluate the backreaction effect for a very simple toy model, namely a Bianchi-I-universe treated as a 
perturbation of a homogeneous and isotropic universe, in section 1101 

The appendix contains a definition of tensor mode decomposition (section[Aj , the equations of motions 
for the first oder complete observables associated to the scalar modes with longitudinal gauge clocks 
(section [B]) and tensor modes (section [C|) , as well as a discussion of issues related to the so called 
linearization instabilities (section ID|) . 



2 Complete observables 

In this section we will give a short introduction to complete observables for finite dimensional systems. 
The generalization to infinite dimensional systems is straightforward. For further details we refer the 
reader to [25l [71 [9] . A very short summary of the necessary details on constraint systems can be found in 
[25] , for a longer introduction see for instance [26] . In order to explain some of the properties of complete 
observables we will express them using a power series, which we are now going to derive. 

To this end we consider a first class (non-degenerate) constraint system with constraints {Cj}JLi. 
Since the system is first class all the constraints generate gauge transformations [55]. Assume that we 
can find a set of phase space functions {T^}^^^ (called clock variables) such that the determinant of 
the matrix 

Af:={T^,C,} (2.1) 

is always non- vanishing (as a function on phase space). In this case we can define an equivalent set of 
first class constraints Ck by multiplying the original constraints with the inverse of the matrix Aj^: 

Ck := C, iA-% (2.2) 

where here and in the following we sum over repeated indices. These new constraints are (weakly, i.e. 
modulo terms proportional to the constraints) conjugated to the clock variables, that is they satisfy 

{T^C4^<5f (2.3) 

where by "~" we denote that this equation holds modulo terms proportional to the constraints. 
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From this property (|2.3p one can proveH that the constraints Ck are weakly Abehan, i.e. their Poisson 
bracket is proportional to terms at least quadratic in the constraints: 

{Ck.Cl} = 0{C^) . (2.4) 

Hence the gauge transformations generated by these constraints commute on the constraint hypersurface. 
This allows us to define gauge invariant observables F[f.j'K-^{T), called complete observables, by a power 
series (assuming that this power series converges): 

F[f;T'<]{r) = E • • {/> CkA, ■■■}, CkA (r^^ - T^O ' ' ' (r^" " T^' ) • (2-5) 

Here / is a phase space function (called the partial observable) and {t^Ij^^i are a set of constants. Due 
to the properties (|2.3l2.4p the Poisson bracket of the series (12.51) with the constraints vanishes at least 
weakly. Hence the complete observable F^.j^k^^t) is (weakly) gauge invariant, that is a (weak) Dirac 
observable. 

Furthermore we have that the complete observable coincides with the partial observable / if restricted 
to the (gauge fixing) hypersurface defined by {T^ — t^}^^i- Therefore we can understand the complete 
observables as gauge invariant extensions of gauge restricted functions, where the gauge is given by 
{T^ = T^}^^i- (In this way one can find a definition of complete observables alternative to the power 
series (|2.5p , see From this observation it follows that associating complete observables to partial 

observables is an algebraic morphism, i.e. we have 



\f-g+h-TK]{''') F[f .rrK]{T) -F^g-rK] + F[,irpK]{T) (2.6) 

as can be also seen in a more elaborate way by using the power series p.Sp . 

The value of the the complete observables is constant on the gauge orbits. In generally covariant 
systems, such as general relativity, the orbits describing time evolution (i.e. foliations of space-times) are 
part of the gauge orbits, or in other words coordinate time evolution is actually a gauge transformation. 
Hence gauge invariant observables are constant in coordinate time. 

Nevertheless one can vary the parameters in the complete observable F[j.y](T). Then the value 
of the phase space function i^[/:T] (''") does change. We can interpret this change as a time evolution 
with respect to the chosen clock variables . Indeed the complete observable F[f.x]{T) gives the value 
of the phase space function / at that "moment" at which the clocks show the values r^. Using 
the representation of the complete observable as the power series p.5|) one can derive the "equations of 
motions" 

d 

which replace the usual time evolution equations. The usual initial conditions are now replaced by 

^[/;T-]M ■ (2.8) 

This closes our short introduction to complete observables. In the remainder of this section we define 
generalizations of energy like observables. 

Here we understand energy as a phase space function which generates the physical time evolution in 
(one of the parameters) r*^. Hence we ask whether there exists gauge invariant phase space functions 
Hm which generate the evolution (|2.7p . that is these functions should satisfy 

^ Fy.T'<]{r)^{Fy,T'<]{r),HM{T)} (2.9) 



Compute {{r^, Cm}, C'i} directly and using the Jacobi identity. Comparing the two results one can conclude that 
KM defined by {Ck.Cm} = fkhl^ 



the structure functions /^m defined by {Ck-,Cm} = fifM^L have to vanish on the constraint hypersurface. 
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where we allowed for a r-dependence of these r-generators. If the variation in say corresponds to 
time evolution (with respect to the physical clock T^) one could interpret the corresponding generator 
Hi{t) as a possibly r-dependent "physical Hamiltonian" [57] (as opposed to the Hamiltonian or scalar 
constraint in general relativity which vanishes on physical configurations). 

Indeed, in the case that the clock variables Poisson commute, it is possible to find functions 
Hk{t) that satisfy ()2.9p . However this equation is not satisfied for all phase space functions /, we have 
to exclude functions which depend on the clocks and functions which do not commute with the clocks 

rpK 

We will look for the gauge invariant functions (t) by assuming that they can be written as complete 
observables F[/j^.yK](T). The partial observables have to satisfy 

d 

— Fy.T'<]{T) ~ ^[{/,Cz,};T'^]W - {F[f:T'<]{T) , ;TA-] (r)} ~ %/,?,^}. ; t^^^t) (2.10) 

where {/,<?}* denotes the Dirac bracket with respect to the gauge {T^ — t^}: 

{/, HlV - {/, hL} - {/, CkUT'', Hl} + {/, T'^UCk, hL} - {/, Cx}{T^, T''}{Cm, Iil} . (2.11) 
For the last equation in (|2.10p we use the general property 

{Fy.T'<] , Flg;T'<] } ^ ^[{f^g}' ; T^] W (2-12) 

which can be proven by using the power series (|2.5p . Here one needs to prove equation (|2.12p only up 
to terms at least linear in the clocks , or in other words on the gauge restricted surface {T^ — t^}. 
Then one can use that the right hand side has to be gauge invariant (because the left hand side is), hence 
it can be written as the gauge invariant extension of a {T^ = r^j-gauge restricted function. 
Comparing equation (|2.10p with equation (|2.1ip we see that the functions have to satisfy 

{/, Cl} {/, hL} - {/, Ck}{T^, Hl} + {/, T''}{CK,hL} - {/, Ca'IIT^, T''}{Cm, h^} + 

0{{T^t)) (2.13) 

where 0{T ~ r) denotes terms vanishing on the gauge fixing hypersurface {T^ — t^}. 

This equation can be satisfied if we chooce = —Pl where Pl is a phase space function such that 
{T^, Pl} — 5^ + 0{C), i.e. Pl has to be a momentum (weakly) conjugated to and to commute with 
the other clocks. Furthermore we have the condition that / has to commute with Pl as well as with the 
clocks {T^}^^^ and that the clocks have to be Abelian. (However the "momenta" {Pl}Y!'=i do not have 
to Poisson commute with each other.) Note that we can add to Hl = —Pl arbitray functions vanishing 
on the constraint hypersurface, for instance Cl- Such additions do not change the associated complete 
observable (on the constraint hypersurface). Also, if one is interested only in the generator associated to 
one specific parameter, say t^, it is sufficient to specifiy Pi. 

Therefore we will assume that the clock variables {T^YJ^^^ are Abelian. In this case one can find 
(locally) symplectic coordinate charts of the form {{T^j^^i, {-Pr-}x=ii {'/'^iPa}} where the Pk are 
conjugated to the and the set {q°,Pa} denotes the remaining symplectic coordinates. 

Then we can define complete observables F[_p^ .j'k-^{t) that generate the evolution in the parameters 

according to 

d 

^F[y.TK](T) ~ {P[^.Tic](T),F[_p^ ;T^-](t)} ^ P[{f,~PLV;T'<]{r) (2.14) 

for / a function of the {q'^,Pa} only. (It is actually sufficient that / commutes with the clocks {T^}^^^ 
and with the momentum H^.) Note however that the complete observables associated to the set {q°',Pa} 
give a complete set of gauge invariant observables. Indeed the restriction to these coordinates as partial 
observables is natural: The complete observable associated to a clock is a constant t^, hence we 
cannot generate evolution in this constant via a Poisson bracket with a generating function. The canonical 
momenta Hk can be solved for by using the constraints. 
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In general it may happen that the Dirac bracket {/, — -Pl}* may depend on the clocks also if / does 
not depend on the clocks. This will occur if the physical Hamiltonian F[_p^ .7-A-](r) is r dependent. This 
is analogous to the situation of an explicit time dependent Hamiltonian in usual classical mechanics. As in 
the latter case we have to add to the time evolution equations a piece that takes care of this r-dependence 
(we will also add a piece, which takes care of the possibly dependence of / on the momenta Pk)- 

d 

.-TR'j (t) ~ {Fy.TK^{T),F[_p^ ;T^-]W} + -^[{/,Pi},T^]W +^[{/,T«}{C'Af,Pi,};T'f]W (2-15) 

holds for arbitrary phase space functions /. 

A T-generator F^-p^ ; t^] (''') does not depend on r if the momentum Poisson commutes with the 
constraints, i.e. if Pl is gauge invariant. This corresponds to a time independent Hamiltonian in classical 
mechanics, i.e. systems with conserved energy. One can always find a set of Abelian clocks {T^}^^^ 
and a set of conjugated momenta {Pk}^^i such that I — min(m, n — m) of the momenta are gauge 
invariant and do not vanish identically on the constraint hypersurface|f| Here 2n is the dimension of the 
phase space. 

We want to emphasize that given a set of clocks the r-generators F^^^p^ . (t) are not uniquely 
determined. To have uniquely determined r-generators (modulo constants and constraints) we have to 
specify the momenta conjugated to the clocks or alternatively the partial observables we want to evolve 
according to (|2.14p . One choice of momenta would be Pl — Cl- In this case the r-generators would 
vanish on the constraint hypersurface. This is consistent with equation p.l4p because functions / allowed 
in this equations have to be gauge invariant, which means that the associated complete observables do 
not depend on the parameters r^. 

In summary we see that complete observables can be used to express dynamics in a gauge invariant 
way and lead furthermore to a generalization of the notions of time (which is specified by using clocks) 
and energy (which is specified by the clocks and a subset of partial observables for which we want to 
describe dynamics). 



3 Approximate complete observables 

We will now develop a perturbation scheme for the complete observables. For such a scheme we need to 
specify which kind of quantities we perturb in, i.e. the quantities assumed to be small. In this approach 
these quantities will be phase space dependent, therefore we will get a good approximation in certain 
regions of phase space. One example of such quantities are deviations from a fixed phase space point 
(serving as a background, e.g. a phase space point describing flat space). We considered a perturbation 
scheme for this case in j9j . 

In this work we are interested in deviations from a whole (symmetry reduced) sector of the phase 
space. In general relativity such a sector could for instance correspond to homogeneous and isotropic 
space-times. In general we will describe the sector, we want to perturb around, by a linear projection 
operator V, which acts on the space of phase space functions. The projection property means that 

v-r^v. 

In particular we can apply the projection operator to some set of symplectic coordinates (which can 
be understood as phase space functions) (x",7ra) with a = 1, . . . ,n where 2n is the dimension of the 
phase space. Then we can write 

^T'-x" + (Id-75) -X" , tt'^ = P • tt" + (Id - P) • tt'^ . (3.1) 

The sector we perturb around is given by the vanishing of the "fluctuations" (Id — P) -x" and (Id — V) -tt". 

^ This statement has to be understood locally, i.e. the phase space function in question may only be defined locally. 
The proof 1281 basically uses that one can locally always find a polarization of the phase space in which (n — m) of the 
symplectic coordinate pairs are gauge invariants |29| . 
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For our purposes we will assume that the projection operator V is of the following forrrQ: There exist 
symplectic coordinate charts (x'^^i i''c)c=i such that the projection on the symplectic coordinates is given 

by 

7^ • x'" = , -P • tt'" = (3.2) 

for a in some subset / of the index set {1, . . . , n} and 

T'-X'^^X'" , P . tt'" = tt'" (3.3) 

for the remaining indices b ^ I. 

The projection operator acts on a general phase space function / by setting all fluctuation variables 
(x'°, 7r'a)ae/ in / to zero, that is 

{V ■ /)(x'", ^'a; x'\ ^'b) - f{V ■ x'\ V ■ 7r'„; V ■ x'\ V ■ 7t\) - /(O, 0; x'\ ^'b) ■ (3.4) 

where the index a takes values in the set / and the index 6 in {1, . . . , n}// 

The fluctuation variables (x'°i t^' a)aei will be considered to be small, that is functions linear in these 
variables are defined to be of first order, functions quadratic in these variables of second order and so on. 
The variables (x'^, T^'b)b<^i are defined to be of zeroth order. 

It is sometimes more convenient not to work with the symplectic coordinates (x' , 7r'c)c=i but with 
some other set (x'^jTTc). The projection operator acts of course also on this set and we can define the 
(degenerate) coordinates 

x^ = {id-r)-x' , Pc = {id-r)-7rc , x'^r-x' , p' = r-7r, . (3.5) 

Since x'^ = + a-nd tTc = Pc + Pc we can expand (suitable) phase space functions in the fiuctuation 
variables (a;°,pa) and introduce a classification of phase space functions by defining the variables {x'^,pa) 
to be of first order. This classification will coincide with the previous one if the coordinates (x''^, 7r'c)c=i 
and (x'^j^'c) are related by a linear symplectic transformation, which we will assume to be the case. 

The coordinates (|3.5p do not need to be symplectic anymore, the Poisson brackets have to be deter- 
mined by 

{X-,Pd}^{V-x',V-n,} , {x^p4-{(Id-^)•X^ (Id-T')-^,} . (3.6) 

The Poisson brackets between fiuctuation variables {x'^,pc) and the "sector" variables {X'^,Pc) vanish. 
Moreover we have the condition on the fluctuation variables that V -x"" = V -pa = 0- Also the coordinates 
{X°-, Pa) will be in general highly degenerate. 

The Poisson bracket of a phase space function of order / with a phase space function of order k can 
in general consist of a term of order {I + k) and of a term of order {I + k — 2) (for l,k > 1). The higher 
order term can arise through the Poisson bracket between the zero order variables {X'^, Pc). The Poisson 
bracket of a phase space function of order I with a zeroth order term is of order I (or vanishes). 

Note that perturbations around a fixed phase space point mo arise as a special case: The projection 
operator is given hy V ■ f — /(toq) • 1, that is V maps all functions to constant functions, where the 
constant is given by the evaluation of the phase space function at the phase space point mg. Symplectic 
coordinates (x''^; c) with properties (|3.2l3.3p can be found starting from any symplectic coordinate chart 
(x^i"''^) and defining x'^ = x'^ ~ x'^("to) and tt^ = tTc — ndmo). The index set / in (|3.2l3.3p coincides 
with the set {!,..., n}. 

For the following we will introduce some notation in order to specify terms of a certain order in a 
phase space function /: with ^'^^ f we will denote all terms which are of order k in /, with I*^!/ we will 
denote all terms in / which are of order less or equal to k. 

''This assumption can be cast into the language of Poisson embeddings, see |16| . It ensures that the kinematics of the 
symmetry reduced system and of the symmetry reduced sector embedded into the full phase space coincide. 
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We define gauge invariant observables of order k as pliase space function wliicli commute with the 
constraints modulo terms of order k (and modulo constraints) . Gauge invariant phase space functions of 
order k can be obtained from phase space functions F which are exactly gauge invariant by omitting all 
terms of order higher than k, i.e. by truncating to I'^IF: 

{WF,Q} = {^^,C,} + {0(fc + l),(")Q + WC7,+...}^0(fc) . (3.7) 

where by 0{l) we denote terms of order I or higher. Here the lowest order term on the right hand side 
will in general appear through the Poisson bracket of the 0{k + 1) term with the first oder term '-^■'Cj of 
the constraint. All other terms are of higher order. 

In particular we can find approximate complete observables of order k by considering their truncation 
to order k. In the following we will assume that the constraints Ck and the clocks can be divided 
into two subsets {{Ch}h€:H, {C/j/gi} and {{T"}h&h,{T'}i€i}, such that the clocks are of zeroth 
order and the clocks are of first order. For the constraints Ch we assume that for an arbitrary first 
order function ^^^J 

{^'^f,CH}^0{l) , (3.8) 

which is satisfied if the constraints Ch do not have a first order term ^^^Ch = 0, however they may have 
a zeroth order term. For the constraints Ci we will assume that the zeroth order terms vanish and that 
the first order terms do not vanish. Consider the power series for the complete observables p.Sp 

oo 

F[/,^.](r) ^ ^_{...{/,(7^J,...,},CKj(r^^-T^^)---(r^'--T^'^) 

oo r 

(r^i - T^i) • • • (t^(--') - T^(-=)) X {t'^ - T^i) • • • {r'' - T^O 

OO OO ^ 

p=0 q=0 ^' 

(r^i -T^i)---(t^'' -T^") (3.9) 

where we used that we can rearrange the constraints in any order and that the clocks commute (at 
least weakly) with the constraints C/. Now set the parameters to zero. Then for the zeroth order 
complete observable associated to a zeroth order function '^"^ we have 

oo 

oo ^ 

- E • • ^"^ChA, ■■■}, (°)C^J(r«^ - T«^) . . . (r^' - T"^) 

q=0 ^' 

(3.10) 

where the second equation holds due to our assumption on the constraints Ch to have vanishing first 
order parts. There only appear zeroth order variables in the second line in p.lOp . hence we can say that 
the zeroth order complete observables associated to a zeroth order functior0 are the complete observables 
of the symmetry reduced sector. The next higher order correction to this complete observable is a second 
order term and can be considered as the correction (backreaction) term to the dynamics of the symmetry 
reduced sector due to deviations from symmetry (in the initial values). 

^ Note that the complete observables associated to a zeroth order function have vanishing first order, hence the zeroth 
complete observable coincides with the first order complete observable in this case. The zeroth order complete observable 
associated to a zeroth order function is therefore a gauge invariant observable of first order. 
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One can also consider for instance the first order complete observable associated to a first order 
function. As we will see these observables describe the propagation of linear perturbations (which are 
linearly gauge invariant) on the symmetry reduced sector. 

Note that this approach allows to find gauge invariant observables to any order k by omitting in the 
series for the complete observables all terms of order higher than k. For this the assumptions we made 
on the clocks {{T^}HGn, {T^}i^i} and the constraints {{CH}HGn, {^i}i&i} ^lot strictly necessary. 
However we will see that with these conditions the computation of the complete observables is similar to 
the usual perturbative calculations involving the "free" propagation of perturbations and their interaction 
as well as the interaction of the zeroth order variables with the perturbations. 

If the power series p.9p for the complete observable converges it defines an exact gauge invariant 
observable which coincides with the approximate Dirac observable ['^li^jj:.^^] (r) modulo terms of order 
{k + 1). If the power series does not converge in some phase space region, this might be due to the fact 
that the clock variables do not provide a good parametrization of the gau^ orbits in this phase space 
region [5S]. In this case one can try to find a set of new clock variables T' , with a better behaviour 
in this respect and such that l|3.8[) is satisfied also for these new clocks. Assume that the complete 
observable Fy,.q^iK^{T'^ ,t'^ ~ 0) associated to these new clock variables and the partial observable 
/' := ['^li^[j.7-jf] (r^, = 0) can be defined. This complete observable will also coincide with ['"'IFj^.t-j (r) 
modulo terms of order fc, as can be seen by examining the power series (j3.9l) for a complete observable 
and using that /' Poisson commutes modulo terms of order k with the constraints. 



4 Application to cosmology 

We will now apply the formalism to general relativity. The sector we perturb around will be the sector 
describing homogeneous and isotropic (FLRW) cosmologies with a minimally coupled scalar field. We will 
work with (complex) connection variables |301 ISlj , however the formalism is independent of the choice 
of variables and can be also applied to real connection variables or the metric (ADM) variables [32] • We 
will follow the conventions in [5T] . 

The canonical variables are fields on a spatial manifold S the coordinates of which we will denote by 
{c°}a=i- We will assume that E is a compact manifold and diffeomorphic to the 3-Torus = S^xS^x 
(in other words the fields are assumed to be periodic). For convenience we will assume that the coordinate 
length of each spatial direction is equal to L = 1. 

The configuration variables are given by a (complex) connection {Ai}j^^i where latin letters from 
the beginning of the alphabet denote spatial indices and from the middle of the alphabet su(2)-algebra 
indices: 

Ai = Ti+(3Ki . (4.1) 

We denote hy (3 = i/2 the Immirzi parameter, is the spin connection for the triads and = 2Kabe'j 
is the extrinsic curvature where is the inverse to the triad e^. The spatial metric can be calculated 
from the triads by qab = ^a^b^jk- The conjugated momenta are constructed out of the triads 

E- = p-'e^^'^'^,n3.« (4-2) 

where e'"*!"^ ejjua totally anti-symmetric tensors with £123 — e^^^ = 1. This gives the following 
relation between the momenta and the 3-metric 

det{q,a)q''' = f3^E;;EiS''^ . (4.3) 
The Poisson brackets between the phase space variables are 

{Aiia),EU<T')}^nSlS'^J{a,a') (4.4) 
where n = SttGn/c^ is the gravitational coupling constant. 
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Furthermore we have a scalar field Lp and its conjugated momentum tt which satisfy the commutation 
relation 

{(^(cr),7r(o-')} = 7'5(o-,cr') (4.5) 

where 7 is the coupling constant for the scalar field. 

The constraints are given by the Gauss constraints Gj ia) , the vector constraints Va (a) and the scalar 
constraints C{a): 

C = K-^p^Fi^,t,uE1,E'^ + ^-\\^:^ + \^}^E'^E]^a^^^,^ + ^3^dei{E^)V{v)) (4.6) 

where V{-) is the potential for the scalar field. Note that we use the scalar constraint with density weight 
two here. In the following it will be convenient to work with the following combination of the constraints: 



Gh 

Da 

S 



- P5iG, (4.7) 
= Va-A^aG, (4.8) 
= G + 2(3''da{E'^Gk5^'') . (4.9) 



In particular Da is now quadratic in the canonical variables and acts as a diffeomorphism constraint. 

We will expand the canonical variables around homogeneous and isotropic field configurations in the 
following way: 

Aa\c7) =Ap5i + aa\o) (3Sl , E-,{<j) = E T^S] + e\{a) p-H] 

ip{cj) =$ + 0(ct) , 7r((7) =n + p{a) . (4.10) 

In this way A, E are real if evaluated on a homogeneous cosmology with flat slicing. This division of the 
phase space into a homogeneous-isotropic sector and an inhomogeneous sector can be implemented by a 
projection operator V acting on the fields defined by: 

Ap^V-Aa' -.^^j^ Aa' da E(3-^=V- E\j ■=\ jji da 

$ ='P-ip If da n =V-TT := / nda (4.11) 



The Poisson brackets between the homogeneous variables and between the fluctuation variables can be 
found by using the projection l|4.1ip : 



{A, i?} - I {aa\<j), e^,(a')} = <-5^-5(a, a') - '^S^S^, 

{$,n} = 7 {^{a),p{a')} =7,5(a,a')-7 ■ (4.12) 



The Poisson bracket between a homogeneous variable and a fluctuation variable vanishes. 

In the following we will raise and lower the indices with the Kronecker symbols S""^ or 5ab respectively 
(and not with the background metric Qat ■— Edat)- 

It will be convenient to work with the Fourier transformed variables, using these the Poisson bracket 
relations simplify. For any field f{a) we define 

f{k) = / expi~ik-a)f{a)da (4.13) 

where k ■ a :— kaa'^ and the wave vector k takes values in 27r The inverse transform is 

/(a)= ^ exp(zA; • a)/(fc) . (4.14) 

fee{27rZ3} 
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Now the homogeneous variables are given by the (|x trace of the) fc = modes of the fields. The Poisson 
brackets for the Fourier modes of the fluctuation variables are 

{0(fc),p(fc')} - 74,-fe' -TWfe'.o (4.15) 

where the additional terms on the right hand side implement that aa°(0) — e°a(0) — 0(0) — p{0) — 0. 

We want to remark that the Fourier transformed variables can be used to define the symplectic 
coordinates used in section [3] in which the projection operator V maps part of the symplectic coordinates 
to zero and leaves the other coordinates invariant. The homogeneous part of the coordinates would be 
given by iV3A, y/3E; $, H). The symplectic pairs that are mapped to zero are given by (aafc(fc), e"''(— fc)) 
and (0(fc), 7r(--fc)) for k ^ Q and 

(\/|(5«ii(0) - «22(0) + ia33(0)) , y|(ie"(0) - 6^2(0) + ^e'^iO))^ 

(^(an(0) - a33(0)) , ^(e"(0) - e'HO))) . (4.16) 

We expand the constraints Cj (where Cq — S, Coa = Da, Cca = Go) in the homogeneous and 
fluctuation variables in order to find the m-th order parts '^'"-'Cj, taking the fluctuation variables as 
first order and the homogeneous variables as zero order quantities. The zeroth order of the constraints 
vanishes except for the scalar constraint: 

(0)5 = K-i6/32^2^2^7-i(in2 + sV($)) . (4.17) 
The first order parts of the constraints are given by 

(n p + E^V'{^) (f) + £;2y($) e^a) (4.18) 

where in the Fourier transformed quantities the partial derivative stands for (daf){k) — ikaf{k). Note 
that the first order of the zero modes of the diffeomorphism constraints Da vanish, this is related to 
the linearization instabilities of backgrounds with compact spatial slices and Killing vectors |34j . In the 
following we will ignore the integrated diffeomorphism constraints Da (0) and show in appendix [D] that 
one can indeed deal with these constraints after one has computed the complete observable with respect 
to all the other constraints. Another way to circumvent the problem of linearization instabilities is to 
couple the system to massless scalar fields and perturb around non-homogeneous field configurations of 
these scalar fields. This will be explained in section [7TT1 

The first order part of the zero mode scalar constraint vanishes too, however the overall scalar con- 
straints has a zero order component. For this reason we do not have a linearization instability correspond- 
ing to the scalar constraint: rather than viewing the second order part of the integrated scalar constraint 
as a restriction on the first order variables, we see it as a correction to the zeroth order part, signifying 
a backreaction effect of the perturbation variables onto the homogeneous variables. This viewpoint is 
possible because the zeroth order variables are part of the phase space, which differs from a perturbative 
approach around a fixed background. 
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Furthermore we need the integrated second order of the scalar constraint: 
('^^(O) = ;^/t-i(2/3^e^'^'^(9„abc(fc) - dbaac{k))e\{-k) + P^E^ {aa^{k)at'{-k) - aa\k)ab%-k)) + 

k 

(3^A^ {e\{k)e\{-k) - e%(fc)e^(-fc)) + 213^ AE {aa''{k)e\{-k) + aa\k)e\{-k))) + 
+ i E^ 5''\dam){db(f){-k) + ($) e\{k) cf,{-k) + 
\ E^V"{^)m<l>{-k) + \ EVi-^) {e\(k)e\{-k) - e%(fc)e^(-fc))) . (4.19) 

To construct complete observables we need to choose a set of clock variables T^{k) where K £ 
{0, Da, Ga; a = 1, 2, 3} and k £27:1? . This has to be done in a way such that at least the zeroth order 
of the matrix 

Af{k,k'):^{T^{k),C,{k')} (4.20) 

is invertible. This will ensure that at least to the lowest order the clock variables provide a good 
parametrization of the gauge orbits (and that is the reason we have to exclude the integrated diffeo- 
morphism constraints, since these start at second order). 

We will choose the clock r''(0) to have a non- vanishing zeroth order part and a vanishing first order 
part. All the other clocks should have vanishing zero order parts. This will ensure that the new constraints 
Cxik) and the clocks have a similar structure as explained above equation (13. 8|) . (We are a bit more 
general here.) The structure of the matrix A is then as follows. The zeroth order of the matrix is of 
diagonal block form, that is 

(°U!^(0,/c') = if J 7^0 or fc' 7^0 

(°U^(fc,0) = if i^T^O or fcT^O . (4.21) 

Since we will quite often need to exclude the index combinations (j — and fc = 0) as well as (K = 
and fc = 0) from the set of indices to sum over, we will introduce indices (j, K, k) to signify that these 
do not assume the values (j = Oandfc = 0) or {K — Oandfc = 0). (Also these indices do not include 
j, K = Da and A: = 0, since we excluded the integrated diffeomorphism constraints.) 

Moreover the first order of the matrix element y^o(0,0) — {T'^(O), Co(0)} vanishes. This structure of 
the matrix A ensures that the constraint 

(70(0) := Y.^,{k!){A-%{k\Q) 

k' 

= Co(0)(^-i)E!(0,0) + Y,C~^{k'){A-%{k',0) (4.22) 

k' 

has vanishing first order. To see this we have to convince ourselves that the first order of (^^^)g(0, 0) 

and the zeroth order of (^^^)g(A;', 0) are vanishing. The latter follows from the fact that A is of diagonal 
block form (14.211) . hence the inverse has the same kind of block form. Furthermore, the first order of the 
inverse of A can be expanded as 

W(^-i)i^(fc',0) = - ^ (°)(^-i)],(fc',fc")«^^(fc",fc"')(°H-4"^)jr(fc'",0) 

k",k"' 

= -Y,^^\A-'yK{k'X)^'^{T^{k"),Com^'\A-y,{m (4.23) 

k" 

For i = and K = Q the sum in the last line collapses to just one term 

«(-4-^)[!(0,0) = -(")(^-i)!](0,0) (i){TO(0),Co(0)} (°)(^-i)[i(0,0) = (4.24) 
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where we used that the first order of Co(0) and r"(0) vanishes. 

Hence the first order of the constraint 6*0(0) vanishes. Also this constraint is the sole one among the 
Cxik) with a non- vanishing zeroth order part: 

^°^CK{k)^J2 '"^^^■(^') ^'\A~%{k',k)^('^Co{0) ^'\A-')l{0,k) (4.25) 

k' 

vanishes for K ^ or k because of the block diagonal form of the zeroth order of A. 

Let us consider the series for the complete observable Fy.xK-^ (r) associated to a function / and with 

parameter values t"(0) — r and (k) — 0: 

F[/^^.](r) ^ ^{•••{/,CK,(fci)},---}^K.(M}(r^Kfci)-r^Hfci))---(r^'-(fc.)-T'''(fcr)) 

00 r 

^ EE E ■(^—^^■■■{f^com^r-s),c^Xk,)},---},c^Xk.)}x 

iT-T%0)Y-' X {-T^^ih))---{-T^^{ks)) 



CXD 00 

EE E 33{---{/>^o(0)},(r-TO(0))«,C^^(fci)},...},C^^(ft,)} 



n\ p. 
P=0 9=0fci,...,fcp 

(_r^i(fc^))...(_r^p(fc^)) (4.26) 

where we denote by {■,-}q iterated Poisson brackets {f,g}q — {{/, 5} and {f,g}o = /. We 
used in the first step that we can arrange the constraints Cxik) in any order since they commute 
up to terms quadratic in the constraints. In the second step we exploited that T'^(O) commutes with 
the constraints Cf.{k) up to terms proportional to the constraints. The result can be interpreted in 
the following way: The complete observable F^f.rpK]^{T) can be calculated by first finding the complete 
observable corresponding to the single constraint Co(0) with parameter value r and then computing the 
complete observable associated to this result with respect to the remaining constraints. One can also 
choose to perform the calculation in the other way around, i.e. first compute the complete observable 
with respect to the constraints Cj^{k) and then to deal with the constraint Cq{0). Here one uses the fact 
that the clocks T^(k) commute with the constraint Co{0) up to terms proportional to the constraints. 



%T-iW - EE E ^{{•••{/'^^K.to---}^*.(M}x 

(-r^^(fci))-..(-T^''(fc,)),C'o(0)| (r-TO(O))' .(4.27) 



q\ p\ 



Assume that / is a zeroth or first order quantity. Expanding the complete observable (I4.26l4.27p up 
to a certain order m in the fluctuation variables, we see that we need the constraints Co(0) up to order 
m for / zeroth and up to order (m + 1) for / first order. We need the remaining constraints up to order 
m for / first order and up to order {m — 1) for / zeroth order . The reason for this is that the lowest 
order in {g, C'^(fc)} is {n — 1) if 5 has lowest order n. However for each Poisson bracket with such a 
constraint the expression gets multiplied with the clock variable T^{k) which is at least of first order. 
On the other hand {(7,(70(0)} is at least of order n if 5 has lowest order n, so one does not loose any 
order in the Poisson brackets with the constraint Co(0). Hence it is crucial that the constraint Co(0) has 
a vanishing first order part, otherwise a perturbational calculation in the usual sense is not possible. 

We can interpret expression (j4.26p in the following way: First we have to evolve the partial observable 
/ with respect to the constraint 6*0(0) which generates time evolution with respect to the clock r"(0). If 
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we want to calculate the complete observable to a certain order m we have to calculate this evolution up 
to terms of order m. This evolution can be broken up into "free evolution" described by the zeroth (for 
the zeroth order variables) and second order (for the first order variables) part of Co{0) and "interaction 
processes" described by the higher order parts (and the second order part for the interaction of the 
homogeneous variables with the inhomogeneities). 

In a second step we have to calculate the gauge invariant extension of the (clock) time evolved function. 
This requires only a finite number of terms in the second (outer) sum since we can omit all terms with 
more than m factors of the (inhomogeneous) clocks (k) . 

In the case that / is of zeroth order, i.e. a "background variable" , the zeroth order complete observable 
corresponds to the complete observable in the isotropic and homogeneous model. Higher (than first) order 
complete observables take into account the backreaction of the fluctuations onto the background. 

For / first order the first order complete observable describes (linearized) propagation of first order 
gauge invariant perturbations. The time evolution of these perturbations is expressed with respect to a 
physical clock (defined by a homogeneous variable). This is different from the usual theory of perturba- 
tions on a fixed background Tf]. 

5 Transformation between different sets of clock variables 

In this section we will explore the dependence of the complete variables on the choice of clock variables. 
Since one can understand the complete observables also as gauge invariant extensions of gauge restricted 
functions this will also enable us to connect different gauges. 

To derive a relation between the complete observables with respect to two different sets of clock 
variables {T^{k)} and {r'^(fc)} we will come back to the interpretation of the complete observables: 
The complete observables F[j.7-jf] (r^) assigns to a phase space point x the value of the function / at 
the point y on the gauge orbit through x at which the clocks (k) coincide with the parameter values 
T^(fc), that is T^ik)iy = T^{k). 

If one replaces in F[f.rpK^{T^) the r^^ (fc) parameters by the complete observables Fjj-a'^/jJ.t/l] (r'^) one 
will get the value of / at that point z on the gauge orbit through x at which the clocks T^{k) coincide 
with the complete observables i^[yK(fc).7-/t] (r' ). Notice that whereas T^{k) changes along the gauge 

orbit the complete observable -Fj^k (r'^) is constant along the gauge orbit, so with the assumption 
that the clocks provide a good parametrization of the gauge orbit the point z is uniquely determined. 
Hence z is a point on the gauge orbit through x which has to satisfy 

T^(fc)|,~F[j,A(fc)^j,,.](T'^) . (5.1) 

The complete observable F\^j.k ^f.yrp,L^ {t'^) on the right hand side gives the value of T^{k) on that point y' 
on the gauge orbit through x on which the clocks T'^ (k) coincide with the parameter values T'^{k). Hence 
the point z in (j5.ip has to coincide with the point y', characterized by the condition T' (fc)|j,/ = r' (fc). 
Therefore we can conclude that 

P[f-a''^]iP[T'<(k)-T''^]{T'^)) - F[f-T'^]i.T'^) ■ (5-2) 

This gives us a relation between the complete observables with respect to two different sets of clock 
variables. If one wants to use this formula in order to obtain F^j.rp,L-^{T'^) from F^j:.'pk-^{t^) one needs 
the functional dependence of the complete observable F^j.^pK] on the parameter values T^{k). However 
with the exception of the parameter t°(0) we set these clock parameters to zero. One can nevertheless 
use formula ()5.2|) if one Taylor expands the left hand side around some fixed values for the parameters 
{k). This would give a formula connecting complete observables with respect to clock variables T'^ {k) 
with complete observables with respect to T^{k). 
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A simpler way to obtain a formula for Fy.j^iL-^ (r' ) as a function of complete observables with respect 
to the clocks T^{k) is to start with the power series for F^j.rp,L^{T''") 

oo 

r=Qki,...,kr 

(r'^^ ih) ^ T'^^ ih)) ■ ■ ■ {r'^^'-iK) - T'^- (fc,)) (5.3) 

where C"A-(fc) = J2k k) and A'f{k, k') = {T'^{k), Cj{k')}. Now we take both sides of 

equation (|5.3p as "input function" / for the complete observable j-kj (r^). The left hand side does not 
change, since the complete observable associated to a Dirac observable is given by the Dirac observable 
itself. Therefore we get 

oo ^ 
r=Oki,...,kr 

(T'^'^ik,) - V-nfcO;T-](r^)) • • • (r'^^-lfc.) - V-'-(fe.);T-](^'')) (5-4) 

as a formula for a complete observable with respect to clocks T'^ as a function of complete observables 
with respect to clocks . 

In a first order approximation (for / first order) and for the case that the clocks r"(0) = r'(0) and 
the corresponding parameter values t°(0) = t'°(0) = r coincide formula (|5.4p reduces to 

^'^%t-](t) ^ ^'^%T-](r)-^(°)F[<o,^^^^,^(^)j^^,,(r)(i)F;^,^(^)^^,,(r) . (5.5) 

k 

Here we set all the other parameter values (k) = t'^ (k) = to zero and we used that r = 
F[T'°(q)-t^]{t)- 

If we have found the first order complete observables associated to some basis of phase space functions 
with respect to one set of clocks T^{k), we can calculate the first order complete observables with respect 
to another set of clocks with the help of formula (|5.5p . Formulae for higher order complete observables 
can be derived by expanding (|5.4p to the appropriate order. 

6 Lapse and shift functions 

As explained earlier the complete observables f[j.Tfc](T^) can also be understood as gauge invariant 
extensions of the phase space function / using the gauge {T^(fc) = r^(fc) VK,k}. Considering the 
complete observables just for one fixed set of parameters would correspond to a "frozen time" picture. 
The time evolution is generated by a constraint, that is time evolution is a gauge transformation. Fixing 
all gauge degrees of freedom would also mean to consider a fixed time. However, we can choose a one- 
parameter family of gauge fixings, as for instance T°(0) = r, T^{k) = for r in (some subset of) R, that 
would represent a varying time. 

A phase space point on the constraint hypersurface gives rise to a solution of the equation of motion, 
that is a space-time manifold. The one-parameter family of gauge fixings defines a foliation of this 
space-time manifold, as well as spatial coordinates on each of the leafs of the foliation. Hence we can 
find (phase space dependent) lapse functions and shift vectors using this foliation and characterize our 
choice for the clock variables and the one-parameter family of r-parameters. 

Lapse and shift can be used to construct the four-dimensional metric using the four-dimensional 
coordinates given by the foliation. This allows one to compare the results of this approach to (covariant) 
methods utilizing gauge fixing. 
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From (|2.7p we see that the (gauge) generator for a translation in the r°(0) parameter is given by 
6*0(0). Hence the lapse functioiO A/"*^ and shift vector J\f^°' can be read off as the coefficients in front of 
the scalar and diffeomorphism constraint. We will also define J\f'^°' as the coefficient in front of the Gauss 
constraint. 

This motivates the definition 

U'{k):={A-%{-k,Q), , (6.1) 

so that we can write 

Co(0) = ^AA^"(fc)Q(-fc) (6.2) 

k 

for our "time evolution" generator (70(0). Note that if one uses a gauge fixing one would restrict the 
consideration to the gauge fixing hypersurface {T^{k) — T^(fc)}, so if one compares (|6.ip to gauge fixing 
one should omit all terms vanishing on this hypersurface. 



7 Clock variables and Hamiltonians 

The clock variables should be chosen such that the zcroth order of the matrix Af {k^ k') is invertible. Once 
one has found such clock variables one can in principle define new clock variables by multiplying the old 
clock variables with the zeroth order of the inverse matrix {A-^)f {k,k'). These new clock variables will 
lead to a new matrix, the zeroth order of which will be given by the identity matrix (with the exception 
of the entry ^o(0, 0) which might differ from 1). 

However these new clock variables might be not very convenient to deal with, since also the first order 
clocks T^(k) will depend on the homogeneous variables, leading to additional terms for the higher orders 
of the matrix Aj^(fc, k') coming from the Poisson bracket between the homogeneous variables in the clock 
variables and the constraints. 

Here we will specify the "inhomogeneous" clock variables T^{k). Hence we allow zero momentum 
k for the Gauss clock T'^'^ but not for the scalar clock T° and not for the diffeomorphism clock T"^" 
(because of the linearization instabilities). One choice for the clock variables is 

^ = £ ^bc 

'j.bciAT„ , LT„ , TL, 



= e'^''%^' etc + etc + 
41 



T" = -W-^dadbe"'' 

= (7.1) 



where W := \J —d'^de- Here we introduced a tensor mode decomposition for the gravitational variables 
(with respect to the flat metric (5ah), the notation is explained in appendix lAl 

This set of clock variables is obviously Abelian and leads to the following zeroth order for the Poisson 
brackets between the clocks and the constraints: 

{r««(fc), (i)Gb(fc')} ^-imih^-w {T^'^k), (i)A(fc')} ^Ed',eb^%k^-k' {r'^"(fc), ^'^S{k')} - 

{r^«(fc), (i)Gb(fc')} - {r^'^(fc), '^'^Dhik')} - E5-,5k^^k' {T'^^k), (i)5(fc')} - 

{r"(fc), (i)G6(fc')} = {r"(fc), ^'^Dtik')} {r"(fc), (^^(fc')} ^-^f^^AE^s^.-k' 

(7.2) 



^ Note that we are using the scalar constraint S here, which has density weight 2. Usually one defines the lapse function 
N-'- as the coefficient in front of the Hamiltonian constraint Cj_ = (-y/g)~^C, which has density weight one. Hence we have 
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By adding to the Gauss clock a term proportional to the difFeomorphism clock we can define new clock 
variables which lead to a diagonal (except for the 5k.-k) factor) matrix fc'). 

^ ^abc(AT+LT+TL)^^^ _ ,abc^^-2Q^Qe TL^^^ ^ W-^Obd" ^^Crfc) (7.3) 

where we assume k ^ 0. (The T'^"^{k = 0) clock is still given by e^^'^Chc-) 
With this new Gauss clock we have 

{r<='^(fc), (i)Gb(fc')} =2l3E5t5k^^k' {T^'^ik), ^'^Db{k')} = {r«"(fc), (^^^(fc')} = 

{T^^ik), (i)Gb(fc')} - {T^'^(fc), ^'^Dbik')} =E5l5k,-k- {T^'^k), (^^(fc')} = 

{TO(fc), (i)Gb(fc')} = {TO(fc), (i)i?,(fc')} = {rO(fc), (i)5(fc')} =-A(3^AEHk,-k' 

(7.4) 

for the zeroth order of the Poisson brackets between the clocks and the constraints. As one can see from 
formula (|5.4p , which gives the relation between complete observables using different sets of clock variables, 
the complete observables using either the clocks (j7.ip or the clocks (|7.3p coincide (for parameter values 
T^(fc) = 0). The reason for this is, that both sets of clocks define the same gauge fixing surface. 

Assume that one has chosen a clock T*'(0), as for instance T°(0) = <&. Then we can define gauge 
invariant functions, that generate the evolution for the complete observables in the r = ''"'^(O) parameter. 
According to section [2] we have to find a momentum Pq conjugated to T'^(O), which has to commute 
(weakly) with the other clocks T^(fc). For the example above we could choose Pq = n. In general we 
will assume that Pq has non-vanishing zeroth order part and vanishing first order part. 

(Since we are only interested in the generator for the evolution in t'^(O), we do not have to specify 
momenta conjugated to the clocks T^(k). However a natural choice would be to choose the first order 
of the constraints '■^•'C'^(fc). For issues arising because of the linearization instabilities see appendix iDl) 

Then we can define the physical Hamiltonian as Ho{t) := Pj/j^ j'K](t) where 

ho^-Po^-Po + Co{0) . (7.5) 
With this physical Hamiltonian we can write 

^F[f,T-]ir) = {Fy^r-]{r),Ho{T)} (7.6) 

for functions / that Poisson commute with Pq and the clocks . If / does not commute with Pq we have 
to add a term ^[{/,Po} , t^^] (''") to the right hand side of equation (j7.6p . This term has the same purpose 
as the additional time derivative dtf that appears in explicit time dependent Hamiltonian systems in 
classical mechanics, where the time evolution equations are given by 

= + . (7.7) 

In general the Hamiltonian Hq{t) will have a non-vanishing zeroth order part (as long as one does not 
choose ^^^Pq = ^^^Cq{Q)) and a vanishing first order part. If one interprets Hq{t) as an energy, this shows 
that also the zeroth order variables contribute to this energy. Hq{t) will be r- independent (i.e. energy is 
conserved) if Pq is a gauge invariant function. This would be the case for Pq = H and constant potential 
for the scalar field. 

Another choice for the clock variables, which is, as we will see in appendix |Bl related to the so called 
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longitudinal gauge pTj is given by 

rj-iGa 



T 



Da 



e { ebc+ etc + ( 



2 



a\de 



1 



e d ■ 



-2/1 T„d 



LLd 



2 

„(1 l^l^r,d \ 

a d— a d) 



(7.8) 



This set of clocks differs from (|7.ip only in the scalar clocks. Notice that the scalar clocks and the diffeo- 
morphism clocks do not commute. The zeroth order parts of the Poisson brackets between constraints 
and clocks are given by 



Jk.-k' 



{T^'^(fc),(i)Gb(fc')} = 
{T0(fc),(i)Gb(fc')} = W-"'5b^fe,-fe' 



{T^'^(fc), (i)i?b(fc')} = E5l5u,-k' 
{T%k),^'^Dbik')} = -W-^AdbSk^-k' 



{T^'^(fc), (i)5(fc')} = 
{T^'^(fc), - 

(7.9) 



For these clocks, which are related to the longitudinal gauge, we cannot give a physical Hamiltonian 
along the lines of section [21 The reason for this is, that the clocks do not Poisson commute with each 
other. However according to equation (|2.1ip . which gives the Dirac bracket, the term that arises because 
of the Non-Abelianess of the clocks is 



k,k' 



{f,CK{k)}{T^{k),T''ik')}{CM{k'),ho} 



(7.10) 



This term is at least of second order if / and ho are gauge invariants to first order. Under these conditions 
a physical Hamiltonian valid for the linearized theory can be defined in the same way as for the Abelian 
clocks (1711). 



7.1 Scalar fields as clocks 

The clock variables we introduced so far are quite non-local, i.e. they require for their definition inverse 
derivatives (in the form of — |fc|~^) or the inverse of the matrix Aj^ requires inverse derivatives. 

One way to avoid this, is to use scalar fields as clocks. Scalar fields are used quite often as clocks, see 
for example [35l [Ml [3 [221 E [37j . As explained in [7j using scalar fields as clocks can lead to huge 
simplifications for the calculation of complete observables. This is related to the fact that scalar fields 
provide a local characterization of spacetime points (as opposed to for instance the longitudinal modes 
of the metric fields, which rather characterize a foliation of spacetime). 

On the one hand we are interested in a homogeneous background, on the other hand we want to use 
the values of four scalar fields to define a coordinate system, i.e. at least the scalar fields defining the 
spatial coordinates have to be non-homogeneous. Hence we will choose these scalar fields to be massless, 
that is having vanishing potential. Since we assumed the spatial manifold S to have topology x x 
we will assume that the scalar fields ip^, A = 1, 2, 3 take values in (parametrized by the values of the 
interval [0, 1[). 

We will use another scalar field tp", taking values in 31, as the time coordinate, that is the clock for 
the scalar constraint. The division into background variables and perturbations is now 



A,'ia)^A[3Si + aa\<j)pSl 
/(a) = $" + (^0(0 , 



■Ko{a) = Ho + po(o-) 
■ka{(t) = + pA{cr) 



(7.11) 
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where ■kq^'ka are the momenta conjugate to the scalar fields Lp^^f"^, respectively. The integrals of (f>^,p^ 
over the spatial manifold are constrained to vanish, this is not the case for the variables (pA, Pa- (For the 
fields ip^^TTA we perturb around a fixed value and not around their averaged value.) 

The Poisson brackets between the gravitational variables are as before (|4.12p . for the matter fields we 
have 

{$",no}=7o, {0°(a),po(a')} = 7o5(cT,a')-7o, {^^{a), pa{ct')} = jAS{cT,a') (7.12) 

where 70, 7_a are the coupling constants for the scalar fields ip'^ , p^ respectively. The matter parts of the 
diffeomorphism and scalar constraints are given by 

matDa = Jo^TTodaip^ + ^ J^^TTAda'P'^ 

A=l,2,3 

= 7o"'no5a(/.° + J2 lA'PA5t + 0{2) 

n^atS = \%\lll+qq'^''dap''d,p° + 2qVo{p°)) + \Y.^~/{l^\+qq'''dap''d,p^) 
^ ^ A 

= + E^V,{^')) + \Y.^a'E^ + 7o"'(noPo + E^V:,[<^°)<P' + i?Vo(<i>°)e\) + 

^ A 

5]7A'(e"W + i?'^''"5a0^) + O(2) . (7.13) 

A 

We could add more scalar fields, these will have the same kind of contribution as the scalar field p'^ . Note 
that the first order of the integrated diffeomorphism constraint does not vanish anymore 

«i?,(A: = 0)= / ^■^'^Da{a)da^ ( ^ lA^RAS^da (7.14) 

SO the problem of the linearization instabilities does not occur for this choice of background (because the 
background values of the scalar fields break the translational symmetry). 

The first order of the integrated scalar constraint is still vanishing, showing that there actually exists 
an exact solution of the equation of motions where all the perturbation variables vanish identically for 
all times. 

Now we can choose as clock variables 

T^''{k)^e"''"'ebc{k) T°''{k) = ^p^{k)5'\ T''{k)^p'^{k) (7.15) 

A 

for all wave vectors = and k ^ Q. As parameter values one has to choose T^°'{k) = 0, T^"-{k) = 
exp{—ika)a°'da and r''(/c) = for fc ^ as well as t'^(O) — t. The zeroth order of the Poisson bracket 
between the Gauss clock and the constraints is as before (|7.2p . for the Poisson brackets between the 
diffeomorphism and scalar clock and the constraints we have 

(°){r^''(fc),A(fc')} = ^b^fc,-fe' , (°){r°(A:),5(fc')} = no4,-fc' (7.16) 

with all the other (zeroth order) Poisson brackets vanishing. Hence the zeroth order of the matrix 
A^{k,k') is invertible on phase space points where llo 7^ and E ^ (from the commutator of the 
Gauss clock with the Gauss constraint). 

These considerations show that one can apply the perturbative formalism also if one uses scalar fields 
as clock variables. Moreover the problem of linearization instabilities does not occur. 
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8 First order perturbations: scalar modes 



Let us consider the first order of the complete observables in more detail. We will assume that we 
are dealing with one scalar field and use the gravitational fields to define the non-homogeneous clock 
variables. Starting from the partial differential equation (j2.7p for the complete observables we will derive 
the equations of motion for the scalar mode perturbations (to first order). The equations of motions for 
the tensor modes are derived in appendix ICl Consider the "time evolution" equation 

-^Flf:T'<]{r) =^[{/,Co(o)};T^]W (8-1) 

where t = t*'(0) is the parameter associated to the clock ^T{k = 0). To simplify the formulae we 
will introduce the notation [[/]](t) :— F[f.xK]{T) and suppress the dependence from the choice of clock 
variables. 

Hence we have 

^ [[{/,Co(0)}]](r) ^ J2 [[{A-')iik',Omr) (8.2) 

k' 

as the differential equation satisfied by the complete observable associated to /. With the introduction 
of lapse and shift fimctions (see section [6|) 

M\k) :^ {A-')ii-k,0) (8.3) 

we obtain the following system of differential equations 

A(i)p(fc)]](^) ^ W[[AA°(0)]](r) W[[p(A:)]](r) + (i)[[AA°(fc)]](r) (°)[[n]](r) (8.4) 

A(i)[[p(fc)]](r) ^ (")[[AA"(0)]](t) (^^[[E^d'^dam) - E'V'{'S>) e\{k) ~ E^V"{'^>) m]]{r) ^ 

(i)[[AAO(fc)]](r) (^^[[E'V'imir) - (i)[[9aAA^"(fc)]](r) (")[[n]](r) (8.5) 

for the first order complete observables associated to the fluctuations <j){k),p{k) in the scalar field and its 
conjugate momentum. Here we assume k ^ 0. 

Now on the right hand side of the equations in (|8.4p there appear also functions of the gravitational 
degrees of freedom, so in principle one would have to add differential equations for these gravitational 
degrees of freedom. However, if we are dealing with only one scalar field in our system, there should be 
only one unconstrained scalar mode degree of freedom. Indeed we can use 

[[C,(fc)]](r)^0 , [[T*(fc)]](r)^0 (8.6) 

to express the (first order complete observables associated to the) lapse and shift functions as well as e^a 
as functions of the (first order complete observables associated to the) scalar field (p and its conjugated 
momentum p and the homogeneous variables. 

For instance for the choice (|7.ip of clock variables we have 

aaT«" = Cabcd^'^^e'' (8.7) 

9„T^" = -W^-2(-|aaa6"e"'' + i5''arf(5,b("e'^'' + ^e'^'')) 

= (-"e"a + iVa) (8.8) 

T° = ^^e^a . (8.9) 



Hence we can use 



,'^\k)]]^0 , [[^e'^\k)]]^0 , [[^^e'^\k)]]^0 . (8.10) 
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Furthermore using the first order of the constraints (|4.18p and the relation (|8.10p gives us 

[[■4^a^^(fc)]]~ + 0(2) , [[^a\{k)]]^^-E-^n^{k) + 0{2) 

7 

[["a%(fc)]] ~ -- (4/32AS2)-i(np(fc) + S'V'($) <j){k)) + -E-^Ilcp{k) + 0(2) . (8.11) 

7 7 

This shows that we can express the gravitational scalar modes as a combination of the matter scalar 
modes. 

The lapse and shift functions (for fc 7^ 0) for the clock variables (|7.ip can be found to be (here we 
abbreviate A/" := ("W°(0)) 

(i)AA^'^(fc) = -UY, *°U~')K"(-fc, k') (iU^(fc', 0) = -AA^ £;-l,5-fc^-fc{T^"(fc'), ^'^5(0)} 

k' k' 

= -NE 2f3^W-^{d'' - i^"^ ^a\) + 0{T) 

(i)AAO(fc) = -AA^(V"i)^(-fc,fc') (iU^(fc',0) = -AA^(-4/32A£;2)-i<5_fc,_fc,{TO(fc'),(')5(0)} 
fc' fe' 
= -M^A-^ +OiT) (8.12) 

where 0(T) denotes terms which vanish with T^{k). Using the equations (|8.1ip the dilferential equations 
for the scalar matter field becomes 

A(i)[[^(fc)]](^) ^ (i)[[AA(p(fc) + i^(Ai?)-itf0(fc))]](r) 

±w[[p{k)]]{T) ^ ^^\[N{E^^^^,m-^fi^'a^m-E'v%^)m + h^i^E)-^'^^pim^ ■ 

(8.13) 

These equations have to be supplemented with the differential equation for the homogeneous variables 
_ff = A, i?, n or <i>, here it is sufficient to consider the zeroth order of this equation: 

^^^\[H]]{T)^^^)[W{H,Sm]] ■ (8.14) 
dr 

Assume that one can find the general solution for these differential equations in dependence on initial 
data for some parameter value t = tq. Since the differential equations (|8.13p is linear in the fluctuation 
fields such a solution for instance for the scalar field can be written as 

Wp(fc)]](r) =.Gi((r-ro);A:;(°)[[i?]](ro)) W[[</.(fc)]](ro) + G,{ir - to); k-^yH]]{ro)) '^'MpikMro) 

(8.15) 

where Gi, G2 can be understood as generalized (free) propagator functions. Now the complete observables 
restricted to the gauge fixing surface have to satisfy 

[[0(fc)]](T)|yo(o)=r,T'i-(fc)=O - '?^(^) 

m]]ir\TOiO)=r,T^ik)=0 ^ H . (8.16) 

For the complete observable associated to the scalar field we therefore have 

(i)[[0(fc)]](T)| Gi ((r - r"(0)); fc; H) 0(fc) + G2 ((r - r"(0)); fc; H) p{k) . (8.17) 

Now we only have to determine '^^^[[0(fc)]](T) away from the hypersurface {T^{k) = 0}. This is done by 
replacing 0(fc) and p{k) in (|8.17p by their first order gauge invariant extensions in T^(fc)~direction 

m -> 0(fc)-^(«)wfc),o^(fc)}T^(fc) 

k 

P{k) - p{k)-J2<^'^{p{k),C,.{k)}T^{k) . (8.18) 
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Hence the first order complete observable associated to the matter scalar mode is given by 
Wp(fc)]](r) G^{iT-T^iO))■,k:H) {m-T.^°Hm,C^ik)}T^ik)) + 

k 

G2{{T-T"{0));k;H) { p{k) -Y,^'^ {p{k),C ^{k)} [k) ) . (8.19) 

fc 

The first order complete observables can be understood to describe the "free" propagation of the 
perturbations on the cosmological background. Higher order complete observables will take care of 
interaction processcfQ given by the higher (than second) order parts of the constraint 6*0(0) as well as 
of backreaction terms arising from the coupling of homogeneous and inhomogeneous variables in the 
constraint. In the next section wc will consider second order complete observables associated to a zeroth 
order phase space function, which capture the lowest order backreaction effects. 



9 Backreaction terms 

Here we want to consider the backreaction effects of the inhomogeneities onto the homogeneous variables. 
To this end we have to find the complete observables associated to a homogeneous variable up to second 
order. (The first order of such a complete observable vanishes.) 

Omitting in (j4.26p with / a function of zeroth order all terms of third order and higher we find 

[2li^[/.y^j(r) ~ (")F[y.T/c](r) + G + / + J (9.1) 

where 



00 

^°^%T-lW = E^{/'^°'^o(0)},(r-r°(0))«=:4;;r(°»(/) and 
9=0 

G - T.{<eT'''\f)^''C^S^^)}{-T^^(k^^)) + 



E ^^°^{{4-e'^"^°"(/)>^'^C'K.(fci)}/'^^/t.(^2)} {'T^Hkm-T^^ik)) ■ (9.2) 

fcl,/C2 

Hence G is the gauge invariant extension to second order with respect to the C^{k) constraints of the 
zeroth order term. The last two terms in (|9.ip are given by 

^ = E E (^H{{//"^Co(0)},„(^)Co},[^lCo(0)},,(r-T°(0))' 

J = ^ Wc^^(fcO}(-T*^(fci)) + 

E '"'{^'H^, ^'^(^^.(fci)}. *'^^K.(fc2)}(-T^^^(fci))(-r^^(fc2)) . (9.3) 

Using the identity 



(gi + 92 + 1)! gi!g2! 



[ ds(i-s)«is«^ (9.4) 
"'0 



^ See 9 for an explicit example of a complete observable taking into account interaction processes in a second order 
approximation around flat space. 
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we can rewrite the term / as 

l-T-T°(0) g 

1= d. ^°^-^^({a}_(/), (^'CoWD (9.5) 

JO 

where we define the free evolution of a higher order term ^"'^ g as 

oo 

«/.ee(^'"^5) - E ^ ^"^3' "^^0(0)}, . (9.6) 

9=0 ^' 

We have also for the free evolution the factorization property a^f^eeif ' d) — '^/ree(/) ' o^/ree(fl')' hence 
the free evolution of a higher order term is determined by the free evolution of its zero and first order 
constituents. Note that for the free evolution of a first order term one can drop all second and higher 
order terms which might arise, i.e. 

(i){(0); (0)^^(0) + (2)6-0(0)} = {(°)/, (°)C'o(0)} (1)5 + ^'^Cm) ■ (9-7) 

In particular the free propagation for the linear terms is linear, i.e. we can express the free evolution of 
a first order term via the propagator functions used in section [S] 
By using that 

{[2lCo(0),r*(i^)} = 0(C)+0(2) and {[2lCo(0), ^^^(fc)} = 0(C2) + 0(2) (9.8) 

we obtain for the term I + J 

rT-T"{o) „ / 

i + J - I d. 4;;r(°)-^)(^ {«}_(/), (2)Co(o)} + 

^'^{Wfreeif), Cq (0)} , « C^^ (fci ) } (-T^H^l )) + 

(9.9) 

The second and third term project out of {aj-^gg(/), (2)C'o(0)} all terms proportional to the clock variables 

T^(k). Furthermore because of l|9.8p one does not need to evolve terms proportional to the hnearized 
constraints ''^^Cj^{k), hence one is left with the evolution of the scalar mode and the two tensor modes. 

In summary we learn that the gauge invariant second order backreaction effect I^li^^y-.j^K] (r) consists of 
two pieces: one is the gauge invariant extension of the homogeneous term F^f-^T''^]{'^) to the appropriate 
order, the other piece comes about through the "interaction" of the homogeneous variables with the second 
order part of the time generating constraint '^^■'C'o(O). Here one needs to consider only the first order 
gauge invariant terms (i.e. the first order physical modes) that arise in this interaction. 

Note that for the second order complete observable I^lFjj.Tic] (r) associated to a zeroth order function 
/ we do not need to consider the integrated vector constraints (related to the linearization instabilities). 
Since these start at second order and t^lF^y.j-K] (r) does not contain a first order term, this second order 
complete observables is already invariant modulo second order terms with respect to the integrated vector 
constraints. 

Higher order complete observables can be calculated by expanding the power series for complete 
observables (|4.26p systematically and by using the identity (|9.4p repeatedly. This will result in a Feynman- 
graph like expansion, that is a sum of terms describing different interaction processes generated by the 
higher order terms of the constraint Co(0). 

In the next section we will calculate for a very simple model the lowest order backreaction effect onto 
the isotropic and homogeneous geometry explicitly. 
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10 Bianchi I as a model for perturbations 



Here, similar to [38] we will consider anisotropic but homogeneous cosmologies of Bianchi-I-type as a 
perturbation of isotropic and homogeneous cosmologies. We will calculate the lowest order effect of the 
anisotropics onto the isotropic and homogeneous variables. The Bianchi I model (with a massless scalar 
field) is solvable, therefore one can compare the results to the exact model. 



10.1 The model 

For the Bianchi I model the connection A^f is given by A,^ — (3 dia,g{Ai, A2, A^) and the canonical 
momentum E'j is given by E'j = /3~^ dia,g{Ei, E2, E3) (see l38j)- Further on we assume that there is a 
massless and homogeneous scalar field $ and its canonical momentum 11. 

The Gauss constraint and the diffeomorphism constraint are trivially fulfilled, the Hamiltonian constraints 
reduces to the following (weight 2 version) 

C = -^{EiAiE2A2 + E1A1E3A3 + E2A2E3A3) + — (10.1) 

K 27 

where (3 — i/2 is the Immirzi parameter and k and 7 are coupling constants. The symplectic structure is 
given by: 

{A,E^}^kS,j {$,n} = 7 (10.2) 

where i, j = 1, 2, 3. 



10.2 Expansion around the isotropic sector 

In order to expand the model around the isotropic sector we define 

E:=^J2^^ (10.3) 

i 

e,:^E,-\Y^Eu . (10.4) 
fc 

This can be implemented through a projector V which projects a phase space function / onto its averaged 
value: 7'(/)(Ai, A2, ^3; i^i, £^2, £^3) = f{A,A,A\E,E,E). We will caU the isotropic variables A and E 
zeroth order variables and the "fluctuations" and first order variables. 
The symplectic structure in these variables reduces to 

{A,E} = ^K, {a„ej} = (5„-i)«:, {<&,n} = 7 . (10.5) 

The fluctuation variables are not completely independent of each other and fulfil the following relations 
(by construction): 

^a. = -0 . (10.6) 

i i 
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In these variables the Hamihonian constraint can be spht into 4 parts: 

C = (3)c+ (10.7) 

(0)^ = ^A^E^ + —n^ (10.8) 

(2) c = M!^£;ya^e^_^^2y g^g_^£;2y„^^^ (^0 9) 

i i i 

(3) c = _M!(AVa,e,e, + £;Va,a,e,) (10.10) 

i i 

("^^C = ^^a.eiajCj . (10.11) 



10.3 Exact solution 



Fortunately this model can be solved exactly and we can compare the exact solution with the one that 
we will derive perturbatively later. We choose lapse N = 1 and can obtain the exact solution by solving 
the following first order system of differential equations. 

A, := {A,C} = 2p^MA,E,+AkEk) (10.12) 

E, := {E,C} = -2P^E,{AjEj+AkEk) (10.13) 

i> := {$,C}=n (10.14) 

ri := {n,C} = (10.15) 

where the indices k on the right hand of these equation are mutually different. (The Einstein sum 
convention does not apply here and throughout this section.) The dot refers to derivative in coordinate 
time t (with the choice iV = 1). One can see that AiEi is a Dirac-observable, i.e. Poisson-commutes with 
C for i = 1, 2, 3. 

We can easily solve this system of differential equations and obtain the following solutions: 

A^{t) = A,exp[-2f3^A,E,t]exp[2(3^^AjEjt] (10.16) 

j 

E,{t) = E,ey.Y>[W^A,E,t\eyiY>[-2(3^^AjEjt] (10.17) 

3 

$(t) = m + $ (10.18) 

n(t) = n (10.19) 

To calculate complete observables we have to specify a clock variable in order to get rid of the physically 
meaningless coordinate time t. For our purposes it is convenient to choose T' = § because it evolves 
linearly in coordinate time: T{t) = t + ^. Inverting this relation and inserting it into (|10.16p . (|10.17p . 
P0.18p . pO.igp leads to the following complete observables: 

F[A.,T^^]{r) = A,exp[-2/32A,i?,(r-|)]exp[2/32^A,i?,(r-|)] (10.20) 

|)]exp[-2/32^A,£;,(r-|)] (10.21) 

(10.22) 
(10.23) 





* ■ 
" n 








^ ■ 
" n - 


|W 


= nr 




^ ■ 
" n - 




= n . 
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For the isotropic variables the observables read as foUows: 

F[A,T=i]{r) = ^ ^ A, eM-2(3'A,E,{r - |)] expp/S^ ^ A,E,{t - |)] (10.24) 

FlE,T=±]{r) = i^£;.exp[2/32A,i?.(r-i)]exp[-2/32^A,£;,(r-|)] . (10.25) 

To be able to compare these exact results with the perturbative calculation later we can Taylor-expand 
(|10.25p around Ei = E, Ai — A. For the complete observable associated to E we obtain: 

$ 

Fie,t=^]{t) = exp[~w(T - -)] X 

x'e+ '^p^a{t - |) E ^^^^ - - E + 

i i 

^f3''u;A{T - ^)2 J2 e^e^ + ^P^uj^{t - ^ a,a, + ^P^cjEir ~ ^ a,e, 

i i i 

+0(3) . (10.26) 

where we introduced the abbreviation uj — AP'^AE and 0(3) denotes terms of order 3 and higher in the 
anisotropic fluctuations. The first order term vanishes due to the condition ^ = ^ = 0. 

i i 

10.4 Using the perturbative approach 

Now that we know the exact solution, we can try to reproduce this results order per order using the 
perturbative approach to complete observables. Here we will only consider the lowest non-trivial order 
of the complete observable associated to the homogeneous variable E. In the first step we have to specify 
a clock function. If we choose T = ^ (in order to fulfil {T, C} = 1) as a clock, the analysis gets as simple 
as possible, because C :— ({T, C})^^C = C. 

In this case the complete observable associated to an arbitrary phase space function / is (formally) given 

by 



= y l^_n^{/, io)c+ Wc}„ . (10.27) 



We can evaluate this sum order per order in the fluctuation variables. For all zeroth order quantities / 
we obtain 

^"^%T=S]W - E J i/' WC}„ =:c.V)^(/) (10.28) 
(i)F[^,^^.j(r) = (10.29) 

. (10.30) 



OO / £N„ 

For first order quantities / we get 

^'^%T.siM = E^^^^'H/'^°'''^^}»-";i(/) (10.31) 

n=0 
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Expression (|10.30p can be interpreted as follows: A zeroth order quantity / propagates with respect to the 
"free" Hamiltonian constraint ^^^C, then there is an "interaction" with the second order Hamiltonian 
constraint ^^^C which yields terms quadratic in the fluctuations. After the interaction zeroth order 
quantities evolve according to the "free" Hamiltonian constraint C and first order quantities according 
to the sum of the zeroth and second order Hamiltonian constraint '^"•'C + ^^•'C (where one can ignore all 
higher than first order quantities that appear in evolving the first order quantity). 

The interpretation of (110. 31|) . which can be seen as the "free" propagation of a first order quantity /, is 
similar. 

We will compute the second order of r^,^ * j (r). Using identity (|9.4p the expression p0.30p can be 
reformulated as follows: 



'■free 



(10.32) 



The first step is to calculate 



a 



(0)C 



a 



(0)C 



(10.33) 



The first two quantities, (|10.33p . are the solutions to the differential equations: 



A := {A, ^°'^C} ^ip'^A^E 
E ■= {E, ^°^C} = -4P^AE^ 



(10.34) 



(10.35) 
(10.36) 



These can easily be found to be 



(A) 



■ (0)C 

*(o)c(£') 

UJ 



Aexp{ujt) 
E exp(— wi) 
AP'^AE 



(10.37) 
(10.38) 
(10.39) 



The second set of quantities can be found by solving 



flfc = {ak,'^'>C}^2p'AEak-2p'A'ek 
ek = {ek,^^^C} = -2(3^AEek + 2f3^E^ak 



(10.40) 
(10.41) 



where the time dependence of the homogeneous variables is given by A{t) 
a*(o)^(£'). The solution to this set of differential equations is given by 



a 



free\0'k) 
t 

free 



(efe) 



UJ A 
= ak exp[ujt] - — (ftfc + —ek)texp[ujt] 
A hj 

u) E 

= efcexp[-tjt] + -(cA; + — afc)texp[-cji] 



q;*(o)^(A) and E{t) 



(10.42) 
(10.43) 
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Figure 1: Bianchi I as a model for perturbations: Here we compare the exact complete observable 
evaluated on a phase space point (Ai, iJ^, $, 11) to the zeroth and second order approximation. The phase 
space point is given hy A = 1, E = Im^, ai = 02 = 0.1, as = —0.2, ei = 62 = 0.1m^,e3 = —Q.2m? 
(where we set the speed of light to c = 1 and the coordinates do not carry units). The momentum 11 
is determined by the constraint (|10.ip . Choosing a value for $ defines the time r for which the above 
values are taken as initial values. 



Now we can calculate (|10.32p step by step: 

1 



Oi 



free 



exp[-cji] 20^A{1 - cut) ^e^e, - 20^ E Ve^a^ + 20^ujt— Va,a 



= -/3^Aexp[-c.(r--)] 



i i 

$,.3^-^ E^ lE^^ . 

" ^2 ^ ^A^ ~ 2 12 "'"'J ^ 

i i i 

1 2 / ^ \ 2 r E sr^ -1 

z i 2 

r 5 ^ — A E ^ — A 3 E"^ ^ — A 



-2w(t - ^) [ 51 + I] + ^ ^ + 



(10.44) 
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This gives the second order contribution 



(2 



n 





(10.45) 



As expected, (|10.45p coincides with (I10.26p . the result we obtained by just Taylor-expanding the exact 
expression. This simple model shows that it is, at least in principle, possible to calculate backreaction 
effects in cosmological applications of general relativity using the methods developed in this work. 

As an illustration we compare in figure[T]the exact complete observable with the zeroth and the second 
order complete observable associated to E. As one can see there the (relative) difference of the second 
order approximation to the zeroth order approximation grows with "time" r, as does the difference of 
the second order approximation to the exact solution. 



11 Discussion and Conclusions 

We presented a gauge invariant canonical scheme for perturbations around symmetry reduced sectors 
of gauge systems. This scheme is applied to general relativiy, in particular to perturbations around the 
cosmological sector. It can be used to calculate the dynamics up to arbitrary high order in the fluctuation 
variables. 

The central objects in this perturbative scheme are complete observables. These complete observables 
are gauge invariant (i.e. Dirac) observables. In the canonical formalism this means that the complete 
observables have to be invariant under (coordinate) time reparametrizations, i.e. constants of motions. 
Nevertheless one can express the dynamics of the theory using the complete observables. To this end one 
has to choose a set of physical clocks. Evolution of dynamical entitities can then be understood as an 
evolution in relation to the clocks. One can even define generalized Hamiltonians, that is gauge invariant 
phase space functions that generate this evolution for the complete observables. 

Different choices of clocks can be interpreted as different setups for a physical measurement. We gave 
explicit formulas relating complete observables associated to different choices of clocks in section [5] In 
particular the clocks define the hypersurfaces over which the averaging (|4.1ip is defined. The complete 
observables evaluated on a phase space point (i.e. evaluated on a spacetime satisfying the Einstein 
equations) give the values of the partial observables, such as the averaged densitized triad E and its 
fluctuations e|J, on the hypersurface determined by {T^{k) = T^(fc)}. Gauge invariance of the complete 
observables ensures that the value of the complete observable does not change if we evaluate the complete 
observable on spacetimes related by diffeomorphisms. From the complete observables associated to all the 
phase space variables (which will give a overcomplete basis of observables) one can calculate the values 
of all kinds of physical observables on the hypersurface {T^{k) — r^(fc)}. For instance one might wish 
to consider instead of the averaged densitized triad some averaged function of the (spacetimqfl) metric 
components: In this case one has first to express the partial observable in question as a function / of the 
zeroth and first order phase space variables and then to consider the complete observable associated to 
this function (expanded to the appropriate order). Note that this complete observable can be expressed 
as the same function / of the complete observables associated to the zeroth and first order complete 
observables. 

*The four metric can be computed by using lapse and shift functions defined in section [B] see also [7]. 
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We considered in more detail first and second order complete observables. Higher complete observables 
can be calculated by similar methods. The zeroth order complete observables coincide with the Dirac 
observables of the symmetry reduced model. The dynamics of the first order complete observables is 
related to the theory of cosmological linear perturbations (see [17] and references therein), the difference 
is that in our case dynamics is expressed with respect to a (global) clock variable, which can however be 
related to choosing a (zeroth order) lapse function. 

Second order complete observables associated to zeroth order partial observables describe backreaction 
effects. We calculated these backreaction effects in section [TU] for a simple model. To allow for an explicit 
expression for these effects in interesting (inflationary) scenarios further approximations, as for instance 
slow roll or long wave length approximations, need to be implemented. This can be done by introducing 
small parameters, which characterize these approximations [39) . The power series for the complete 
observables has then to be expanded to a certain order in the fluctuation variables and to a certain order 
in these new parameters. 

A perturbative scheme around a whole phase space sector has several advantages compared to an 
expansion around a fixed phase space point. Firstly, some of the degrees of freedom, namely these 
corresponding to the symmetry reduced sector are treated non-perturbatively. These degrees of freedom 
are used to define the clock in relation to which we express the dynamics of the theory. Compared to 
defining a field theory on a fixed background, where the background is used to define time, the clock in 
this approach is a fully dynamical entity, as one would expect in the full theory. Indeed one would expect 
problems with gauge consistency to higher order if one uses a background time, since such a time cannot 
take into account backreaction effects onto the physical clock. 

Also a quantization of the theory should not only quantize the perturbations on a fixed background 
but rather consider the quantized perturbations on a quantized geometry, described by the (quantized) 
symmetry reduced sector. The variables describing the symmetry reduced sector contribute for instance 
to the Hamiltonian generating the physical time evolution. From this perspective it might be fruitful to 
reconsider certain concepts, like energy and vacua, from quantum field theory on curved spacetime. 

As in the case for perturbations around a fixed background, where the background is given in certain 
coordinates, the symmetry reduced sector is described by using a certain type of coordinates. That is, 
there exist configurations which are not included in the symmetry reduced sector, but which are never- 
theless physically symmetric. These configurations can be obtained by a "non-symmetric" coordinate 
transformation from a symmetric configuration. Degrees of freedom gained in this manner correspond to 
gauge degrees of freedom. 

However by expanding around a sector describing a family of solutions in different coordinate sys- 
tems one can investigate the non-perturbative effects of such gauge degrees of freedom compared to an 
expansion around a smaller sector (which could be just one phase space point), which would treat some 
of these gauge degrees of freedom perturbatively ^40j . That is we can embed symmetry reduced models 
into each other and explore in this way the reliability of symmetry reduction. 



Appendix 

A Tensor mode decomposition 

Similar to the longitudinal and transversal modes for a vector field on one can introduce tensor modes 
for a tensor field. For a proof of the completeness of these modes, see [33]. To begin with we define the 
projector onto the transversal modes of a vector field by 

{p-v)a:^pl-vt,:^{5'', + W-^ ■dadb)-vt, . (A.l) 
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This allows us to introduce the following basis of tensor modes: 

(^^P • T)ab = {6a - pV) ■ pt ■ Ted 2 left long, right transv. modes 

(^^P • T)ab = [51 - pI) ■ {5i ~ pt) ■ Tad 1 left and right long, mode 

C^P ■ T)ab = pI ■ {St - pt) ■ Ted 2 left transv. right long, modes 

p . T)ab = \pab ■ p'^'^ ■ Ted 1 symm. transv. trace part mode 

lATp , y^^^ ^ 1 . „ , , j,^^ ^ antisymm. transv. mode 

{STTp . 2^)^^ ^ 1 (-pC .pdj^pc.pd^ PabP"'^) ■ Ted 2 symm. transv. tracefree modes (A.2) 

Using the projector property p • p = p, it is easy to see that the projectors ^ P are orthogonal to each 
other and satisfy ^ P ■ ^ P = S-^^ ^ P. Furthermore the set of projectors is complete, that is 



X 



■PaS^m • (A.3) 

We will denote the tensor modes by ''^Tab '■— P ■ T)ab- 

B First order perturbations: scalar modes in longitudinal gauge 

Here we will derive the equations of motions for the first order scalar perturbations with the choice (|7.8p 

pGa ^ ^abc^^^ 

= e^'^i^^Cbe + '^^Cbe+'^'^ebe) 

T^" = -W-\-^W-^d''ddde + ^d^5de-de5^-ddS:)e''^ 

= -W-^i-d"" ^^e^'d + ^<9" ^e'^d - de ^^e""^ - dd ^^e'*'^) 

= W-''{\5^'^ + lW-^d^d'^)aed 

= W-W^d^d-'^'^a'^d) . (B.l) 

for the clock variables. As we will see this choice is related to the longitudinal gauge: In the notation of 
section [8] we have 

\\T^-]\ = W^W-^i-lW-^d'^ddde + \d''bde - de^l " ddSDe^^W ^ . (B.2) 

(Note that the brackets [[ • ]] now refer to complete observables with respect to the clock variables (|B.1[) .) 
We want to translate this condition on the triad perturbation to a condition on the spatial metric qab- 
With the relation 

det(g)g'''' = ii^E^^E] (B.3) 

between the densitized triad variables E^ and the inverse metric q"-^ (where Aet^q) is the determinant of 
the metric) we get 

qab = E^ab + e^eSab " ^ab " ^ba + 0(2) =: E^ab + Kb + 0(2) . (B.4) 

Hence Cab + e.ba = —hab + hedS'^'^Sab + 0{2) and together with [[cab]] — [[eba]] from the Gauss clock (in 
IB.ip we obtain 

2 (i)[[r^-]] = (i)[[Ty-2(-iT^-2a-5rf9e + ^d'^Sde - dJ2 - ddS^e){h''^ - S^^K)]] 

= (l)[[M/-2(_lM/"25-9^a^ + ^d'^Sde - deS2 " 9d<5,")/i''l] 

~ . (B.5) 
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This is the longitudinal gauge condition for the three metric (see [17]). Furthermore we can calculate the 
lapse and shift functions introduced in section [6l 

(i)[[AA^"(fc)]] ~ 

^^^iJ^'^'^ik)]] ~ (i)[[AA5V^,(fc)]] (B.6) 

where we used the definition A/" = '•"W°{0) = ((°){r"(0), C'o(0)})-i. In particular we see that the shift 
function [[A/'^°(/i;)]] vanishes, which is the other gauge condition in the longitudinal gauge (without 
coupling to spin one (vector) matter fields). Hence we can see our first order observables as gauge 
invariant extensions of longitudinal gauge restricted functions. 

These lapse and shift functions can be used to derive in the same way as in section [8] the equations of 
motions for the scalar matter field *^^^[[0]](t): 

A(i)[[p]](^) ^ (^)[[N-{E^d^da(P~^E^V'{^)e\-E^V"i<^)(b+lE^V'i^)e\mT) . 

These equations lead to the following wave equatior[^ for the scalar mode: 
^Wp]](r) c (")[[AA]](r) A(o)[AA]](r)A(i)p]](^) + 

(i)[[AA2(s2^°a,0- ii;V($)e"a - £;^r'($)0)]] - 

(B., 



(B.7) 



Using ^^'>[[T^ik)]] = + 0(2) and [[C'^(fc)]] ^0 + 0(2) one can replace the metric mode ^^'>[[e\]] 
in (jB.Sp by some combination of the scalar field '^■'[[0]] and its first r-derivative, however this does not 
lead to a simple equation. 

Another possibility is to derive a wave equation for the metric mode ^^-'[[e'^a]] in the same way as for 
the scalar field mode: 

A(i)[[e-J](r) ^ (^^[[-mf3'E'a\]]{T) 

A (i)[[a\]](r) ~ ^'^[[MiS-6ae% - 4/?^^ e^, + 8(3^AEa% + 3^i?V(a>)0)]](r) 



(B.9) 



Hence we obtain the wave equation 



^W[[eM](r) (")[[AA]](r)^(°)[An](r)^W[[e^]](r) + 

^^^[[~AfH/3^E^{d''dae^c~'iP^A^e^, + 3^E^V'{^)(f>)]]{T) (B.IO) 
where we can replace the matter scalar field by 

(i)[M](r)^(i)eH-iiAe^,]](r) + (°)[[|(AA4/32£;)-i]](r) . (b.H) 



9 For conformal time , that is if one chooses the clock T°{0) in such a way that ^"W = [{T0(0), °)Co(0)}]-i = E'^ 
the wave equation IIB.8I 1 coincides with the wave equation for the (first order gauge invariant) scalar field in 1171 . (The 
metric scalar mode ^p{k) used in il7l can be computed to be i/) = —^E~^e'^c- Also the (first order) gauge invariant Bardeen 
potential gl] is given by ^ = '-^^[[-^E-^e'^c]]-) 
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Using the time evolution of the homogeneous variables A, E, (f>, 11 and specializing to conformal time 
{J\f = E^^) we therefore have as the wave equatiorP^ for [[E'^^e'^a]] 

^W[[E-'e\]]{r) ^ W[[4f3'd^daE-'e%]]ir) + -^('^^[[E]]ir)^'^[[E~'e%]]{r) + 

(«)[[2n-i]](r) A (0)[[n]](.) ^ W[[E-'e^Mr) - 

(°)[[i?-in-i]](r) A(o)[[i;]](^) d (o)jp]](^) W[[E-^e%]]ir) . (B.12) 



C First order perturbations: tensor modes 

Here we will derive the differential equation for the tensor modes of the gravitational field, that is the 
STT-modes. The STT modes ^'^'^a"'^ and ■S^^'gob already gauge invariants of first order, since they 
commute with the first order part of the constraints. Hence according to equation (|5.5p which connects 
the first order complete observables with respect to different choices of clock variables, the first order 
complete observable associated to the STT-modes are independent from the choice of clock variables 
T*(fc). 

Indeed in the differential equation (|8.2|) for the first order complete observables associated to the STT 
modes (where A^-^ = ^°\T°{0), S{Q)}) 

- (^)[[{^^^a.b, C'o(0)}]](r) ^ (i)[[AA{^^^a,,, S(0)}]](r) 

- (^)[[{^^^e.,, Co(0)}]](r) ^ (i)[[AA{^^^e,,, S(0)}]](r) (C.l) 

all terms which may depend on the choice of the T^(fc) variables drop out. Here we used the definition 
(|4.22p of Cq (0) and the fact that the 5TT-modes commute with the first order part of the constraints. 
With the second order scalar constraint (|4.19p we get 

^%N{2pEDil s^^af, - 2{P'A' + '^E^Vm ^^^e,, + 2P^AE ^^^a,,)]] 
^ (i)[[AA( - 2pEDil ^^^e/e + 2(3'E' ^^^a,, - 2^^ AE^^^e,,)]] (C.2) 

where D^J^ = ^e'^'^'^dc{5dbSci + ^da^l)- If D acts on S'TT-modes it simplifies to D^^l — eb'^'^dcSl, moreover 
on S-TT-modes we have {D ■ D){l = -d^dc SfJl 

We want to derive a wave equation for the ^"^^Cab niodes. To this end we have to calculate the second 
(r-) time derivative of ^^^Cab- In the process we have to take also into account the r dependence of the 
homogeneous variables and to replace the '^"^"^aab-modes by the r-derivative of the ^"^"^Cab niodes with 
the help of the second of the equations (jC.2p . The resulting equation of motion for the ^'^'^Cab modes is 

(!)[[_ AA2(4/32^2g=^^ + 4/32£;3£y($) + l&p^A^E'^f'^'^ Cab]] ■ (C.3) 

The wave equation simplifies if we consider instead of ^^^Cab the rescaled variable E~^ ^"^^eab- (The 
correspondence to the metric variables is given by E~^^'^hab = —2E~^^^^eab where hab is the deviation 

This wave equation coincides with the wave equation for the Bardeen potential ^ = [[i?~^e'^c]] in |17| . 
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of the spatial metric qab from the isotropic and homogeneous background qab = E6ab + hat ='■ Qab + hab- 
Hence the rescaled variable corresponds to {—2)Q"^'^'^hcb-) 

(C.4) 

Specializing to conformal time (J\f = E^^) wc find 

^^^eab]] ^ -n[E-']]-^nm^ ^'\[E-' '^^eab]] - 4l3^d^d,(%E-' ^^^e^b]] • (C.5) 



D Linearization Instabilities 

So far wc left out the discussion of the integrated diffcomorphism constraints, which start with second 
order terms (if one does not use three massless scalar fields as clocks for the diffcomorphism constraints). 

Here we will show for a general first class constraint system, that one can apply the complete observable 
method with respect to a siibsct of the constraints, also if these do not form a subalgcbra. In our case 
this subset is given by all the constraints with the exception of the integrated diffcomorphism constraints. 
After constructing these (partially invariant) complete observables one can use these to get fully invariant 
observablcs. 

We will explain the procedure for a finite dimensional system, the generalization to field systems is 
straightforward. To start consider a first class constraint system with constraints {(7j }^j and subdivide 
this set into two subsets of constraints {{Ca}^=i, {C'a}^=i^}- We want to apply the complete observable 
method to the first subset {Ca}^i of the constraint set. Hence we choose a set of clock variables {T^}'^^-^^ 
and define in the same way as we would deal with the full set of constraints the new constraints 

CA-.^iA-'TACa where A^ = {T^,C„} (D.l) 

These new constraints satisfy 

{TB,CA} = 5'i+\TCc (D.2) 

for some phase space functions A^*^. 

Furthermore we want to add to the constraints {Ca\ another set of constraints such that we get a 
complete system of constraints and such that these added constraints commute at least on the constraint 
hypersurface with the clock variables {T^}. Therefore we define the constraints 

Co, ■= a-{T^,C„}C'^ (D.3) 

which satisfy 

{T^,a} = 0{C) . (D.4) 

We want to consider the Poisson algebra of the new constraints Ca and Ca- To this end we define 
the structure functions / by 

{Ca.Cb} = 

{CA:Cf}} = 
{C^,C/3} = 

Then using the Jacobi identity we can calculate 

0{C) = {Ca,{Cb,T^}} = {T^ACb,Ca}} + {Cb,{Ca,T^]} 

= {T^, JEaCd + f^AC'a} + {Cb, -5^ + 0{C)} 

= fBA + 0{C) . (D.6) 



/aB^C + fAB^a 

f^aCc + taCa ■ (D.5) 
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Hence we have f^g — 0{C) and therefore 

{CA,CB}^0{C^) + f2BCc . (D.7) 

Similarly 

{T^,{Cb,C^}} = {T^JEcCc + fLCp} 

= fL + 0{C) . (D.8) 

On the other hand 

{T^,{Cb,C^}} = {CB,{T^,C^}} + {Ca,{CB,T^}}^OiC) (D.9) 
hence fg^ = 0{C) and we have that 

{Cb,C^} = 0{C^)+~fl^C^ . (D.IO) 
In the same way one can proof that 

{C^.Cp} = 0{C^)+~flpC, . (D.ll) 

In summary we learn that the constraints {Ca} are weakly Abelian modulo terms proportional to 
the constraints Ca and that the constraint set {Co] forms an ideal (modulo terms quadratic in the 
constraints) of the full constraint algebra. 

Assume that we have a phase space function / that is (weakly) invariant under the constraints {Ca\- 
Then we can use the power series for complete observables p.Sp just with the constraints {Ca] and 
compute the complete observable associated to /: 

oo 

F[f,T^^{r^) =^ E fc! {• • • {/' J, • • • }, J(r^^ - T^^) • • • {r^^ - T^") . (D.12) 

k=0 

Because of the properties (jD.7ID.10|) the resulting function is (weakly) invariant under both sets of 
constraints {Ca} and {Ca}- (Just consider the Poisson bracket of (jD.12p with a constraint Ca and Ca-) 
That is to find a fully gauge invariant observable we have to start with a partial observable / that 
is invariant under the constraints {Ca}- However it is also possible to calculate complete observables 
(jD.12p associated to non-invariant functions / and then to find fully gauge invariant observables using 
that 

F[h-h+h;T-]ir^) ^ %;T-](r^)-%;T-](r^)+%;T-](r^) . (D.13) 

Hence we can obtain a fully gauge invariant observable by combining complete observables associated to 
non-invariant partial observables fi algebraically such that the same algebraic combination of the fi is 
invariant under the constraints Ca- 

Applied to our situation this means that we can ignore the integrated diffeomorphism constraints in 
computing complete observables. However in the end we have to find combinations of partial observables 
that are invariant under the (altered) integrated diff'eomorphism constraints CoaiO) defined by 

CoaiO) = CoaiO) - Y,{T''{k), Coam CK^k) . (D.14) 

(The sum does not include the values {K, k) — {Db,0)-) Note that also the altered integrated diffeomor- 
phism constraints start at second order. 

Since the integrated diffeomorphism constraints generated coordinate translations, functions of the 
form f{k)g{—k) for some k are invariant under the constraints CoaiO)- Here we denote by f{k) the 
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Fourier transformation of a field variable /(cr). Therefore such functions are a good point to start to look 
for functions which are invariant under the constraints (jD.14|) . 

Concerning the discussion of the gauge invariant r-generators Ho{t) = F[_pg j^k^{t) in section [71 
the results do not change if we demand that / in equation (j7.6p and Pq arc invariant (to the order in 
question) under the altered integrated diffeomorphism constraints CoaiO)- 
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